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ABSTRACT 



A digital computer program for the geometrically nonlin- 
ear static and dynamic response of arbitrarily loaded shells 
of revolution known as SATANS is modified, removing a re- 
quirement that the applied loads be symmetric about a datum 
meridional plane. The ability to include arbitrary 
naturally-occurring initial geometric imperfections of any 
shape, location or extent is added to the program. Two 
specific axially loaded circular cylinders are analyzed 
using actual initial imperfections. Detailed numerical 
studies are conducted to determine the buckling loads of the 
perfect and imperfect cylinders. Several sets of boundary 
conditions are considered. An analysis procedure for imper- 
fect shells is proposed. The results indicate that a means 
to conduct a numerical analysis of an arbitrarily imperfect 
shell of revolution has been developed. 

A study of the buckling behavior of the clamped shallow 
spherical cap modelled under various boundary conditions at 
the pole is included as an addendum. 
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TABLE OF SYMBOLS 



a 


= reference length 


b 


= nondimensional inplane stiffness 


E 


= Young's modulus 


E o 


= reference Young's modulus 


h 


= thickness 




= reference thickness 


H 


= rise of a shallow spherical cap 


m 


= harmonic number, meridional 


M s' V M s0 


= bending and twisting moments per unit 
length 


m s • m 0' m S0 


= nondimensional bending and twisting 
moment series coefficients 


n 


= harmonic number, circumferential 


N s' N 0' N S0 


= membrane forces per unit length 


p 


= nondimensional pressure load series 
coefficient 




= effective transverse force 


q> q s » qe 


= normal, meridional and circumferential 
components of applied pressure loading 


q/q 0 


= ratio of applied load to critical load 


$0 

w 

•>* 

CD 


= principal radii of curvature 


r 


= normal distance from the axis of the shell 


s 


= meridional shell coordinate 


V fc 0' t s0 


= nondimensional membrane force series 


t T 


= nondimensional thermal membrane force 
series coefficient 


U, V 


= displacements tangent to the meridian and 
tangent to the parallel circle 
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U, V 


= nondimensional series coefficients of U, V 


W 


= displacement normal to the reference 
surface 


W 


= displacement normal to the reference 
surface due to initial imperfections 


w 


= nondimensional series coefficient of W 


8/ 8 S , 8 e , 8 sQ 


= nondimensional nonlinear term series 
coefficients in the strain-displacement 
equations 


A 


= distance between stations 


e s' e 0 ' G s0 


= nondimensional reference surface strain 
coefficients 


£ s' e 6' £ S0 


= nondimensional reference surface imperfec- 
tion strain series coefficients 


S 


= shell coordinate perpendicular to the 
reference surface 


0 


= shell circumferential coordinate 


n, n s , Hq, n s g 


= nondimensional nonlinear term series 

coefficients in the equilibrium equations 


X 


= nondimensional geometric descriptive 
parameter for a spherical cap 


V 


= Poisson's ratio 


a o 


= reference stress level 


p 


= nondimensional radius 


u s' w 0 


= nondimensional curvatures 


*' $ S' $ 0 


= reference surface rotations 


V ? 0 


= reference surface rotations due to initial 
imperfections 


<K <t> S ' 4>e 


= nondimensional series coefficients for 
reference surface rotations 


4>s ' $6 


= nondimensional series coefficients for 

reference surface rotations due to initial 
imperfections 


/"q_7, /“a 7 


=4x4 boundary condition matrices 
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{$,} =1x4 boundary condition matrix 

/ — E_7f / F_/, / G_/ = matrices / Refs. 3 and £/ 



£ = a summation, on the harmonic index num- 

i ber n, where n = 1, 2, •••, 00 



SUPERSCRIPTS 

1 = indicates the coefficient of a previously 

existing term 

2 = indicates the coefficient of a newly intro- 

duced term 

12 = indicates the coefficient of the product 

of a previously existing term and a newly 
introduced term 

(n) = Fourier index 

" = (prime) , partial differentiation in the 

meridional direction 

• = (dot) , partial differentiation in the cir- 

cumferential direction 
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I. INTRODUCTION 



Shells are used as structural elements in nearly every 
space vehicle, aircraft or other structure where light 
weight is of importance. Consequently, a large number of 
digital computer programs for the structural analysis of 
this shells have been developed. Many of these programs 
were briefly described by Bushnell in Ref. 1. An assessment 
concerning the analysis capability based on computational 
systems that have application to a broad class of shell 
structures was presented by Hartung in Ref. 2. In his as- 
sessment, Hartung discusses the limitations, scope, accuracy 
and efficiency of virtually every major shell structural 
analysis program currently available and under development. 
Of particular interest here, is the computer program SATANS 
— £tatic And Transient Analysis, Nonlinear, £hells — de- 
veloped by Ball for the geometrically nonlinear analysis of 
arbitrarily loaded shells of revolution /Refs. 3 and Aj . 

First, in section II, a brief description of SATANS is 
presented including the basic theory, the program capabili- 
ties and limitations, the solution method and procedure and 
some general comments. 

Second, in section III, SATANS-I is constructed. Modi- 
fications to SATANS making the program more user oriented 
are discussed. Included in these modifications are a reor- 
ganization of most of SATANS into a 'storable package,' 
incorporation of an equipment-independent plotting package 
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which uses a high-speed line printer, modification of itera- 
tion control, discussion of extrapolation technique and 
alteration of boundary load application procedures. 

Third, and most important, the subject of imperfections 
is treated in section IV. Thin isotropic shells of revolu- 
tion have consistently displayed experimental buckling loads 
much lower than the theoretically predicted buckling loads. 
This disagreement has generally been attributed to (1) , the 
mathematical inability to correctly model the boundary 
conditions as they exist on the test specimen, and (2) , a 
high degree of imperfection sensitivity, i.e., naturally- 
occurring imperfections in the initial shape of the shell 
cause a premature triggering of critical buckling modes and 
early shell collapse. An attempt by Dantone /Ref. 5 / to 
modify SATANS and to include initial imperfections in the 
analysis of cylindrical shells was severely hampered by a 
restriction on the use of SATANS, namely, all applied 
loading and initial conditions — and hence, initial imper- 
fections — must be symmetric about some meridional plane. 
Discussion of the removal of this restriction and introduc- 
tion of initial imperfections into the program are covered 
in this section. 

Fourth, and also in section IV, the implementation of 
the initial imperfection capability is presented. The 
resulting program is called SATANS-II. Steps for actual 
construction of SATANS-II from SATANS-I are presented as 
are instructions for its use. 
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Fifth, in section V, actual imperfect shell analyses 
are performed on two specific axially loaded circular cylin- 
ders. The actual topography of the cylinders has been 
measured by Arbocz and Babcock /Ref. 6 / and this data pro- 
vides the opportunity to include actual naturally-occurring 
initial imperfections into the analysis by SATANS-II. 
Detailed numerical studies are performed to determine the 
classical and consistent bifurcation buckling loads of the 
perfect cylinders and the results are compared with previ- 
ously published results. Stability solutions for the 
imperfect cylinders under six boundary conditions are 
obtained and compared with experimental values. An analysis 
procedure for imperfect shells is proposed based upon these 
results . 

Lastly, as an addendum to this report, the clamped 
shallow spherical cap is investigated under various boundary 
conditions at the pole. While on an industrial tour with 
the Structural Analysis Division of the Palo Alto Research 
Laboratory of the Lockheed Missiles and Space Company, 
Sunnyvale, California, Ryan noted that few shell structural 
analysis programs have the capability of analyzing a shell 
of revolution which has an initial or a final pole — 
continuous structure through the axis of revolution. This 
investigation attempted to determine the effects on the 
solution caused by substitution of a small hole or a small 
rigid plug for the pole. 
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II. DESCRIPTION OF SATANS 



The computer program known as SATANS is capable of per- 
forming static and transient analyses of thin, linearly- 
elastic shells of revolution subjected to arbitrary loads 
and temperature distributions. It is based on Sanders' 
nonlinear field equations /Ref. ~fj for the conditions of 
small strains and moderately small rotations. The program 
can be used under the following conditions: 

1. The geometric and material properties of the shell 
must be axisymmetric , but may vary along the shell meridian. 

2. The shell material must be isotropic, but the 
modulus of elasticity may vary through the thickness. 

3. Poisson's ratio must be constant. 

4. The boundaries of the shell may be closed, fixed, 
elastically restrained or free, and loaded or unloaded. 

5. The applied forces, pressures, temperatures and 
initial conditions must be symmetric about a datum meridional 
plane . 

A. GEOMETRY AND COORDINATE SYSTEM 

The coordinate system used to locate points within the 
shell is the s, 0, c; orthogonal system, where s is the 
meridional distance along a reference surface, 0 is the 
circumferential angle measured about the axis of revolution 
and £ is the perpendicular distance from the reference 
surface. The shell geometry and coordinate system are shown 
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in Figure 1. Positive directions for displacements, rota- 
tions, forces, moments and loads are shown in Figure 2. 

B. METHOD OF SOLUTION 

The governing partial differential equations are reduced 
to an infinite number of sets of four second order differen- 
tial equations in the meridional and time coordinates by 
expanding all dependent variables in a sine or cosine series 
in terms of n0 , where n is the Fourier index (harmonic) and 
0 is the circumferential coordinate. The sets are uncoupled 
by utilizing appropriate trigonometric identities and by 
treating the nonlinear coupling terms as 'pseudo loads.' 

The meridional derivatives are replaced by conventional 
central finite difference approximations, and the displace- 
ment accelerations are approximated by the implicit Houbolt 
backward differencing scheme /Ref. 8/. This leads to sets 
of algebraic equations in terms of the dependent variables 
and the Fourier index or harmonic. At each load or time 
step, an estimate of the solution is obtained by extrapola- 
tion from the solutions at previous load or time steps. The 
sets of algebraic equations are repeatedly solved using 
Potters' form of Gaussian elimination /Ref. 9/. The pseudo 
loads are recomputed after each solution and the procedure 
continues until the solution converges. 

In the static analysis, a load-displacement history is 
obtained by an automatic incrementing of applied loadings 
in equal steps until the number of iterations required to 
achieve convergence in any harmonic exceeds a prescribed 
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SHELL GEOMETRY AND COORDINATE SYSTEM 

FIGURE l 
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A. DISPLACEMENTS AND ROTATIONS 




POSITIVE DIRECTIONS FOR VARIABLES 
FIGURE 2 
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limit. The load increment is then reduced by a factor of 
five and the procedure continues until a prescribed number 
of load steps, or a prescribed number of load step reduc- 
tions have been made. 

Since the method of solution is based on the pseudo load 
approach, the shell reacts equally, in a linear fashion, to 
any change in either the applied load or the pseudo load. 
Thus, failure of the solution to converge in any harmonic 
can be attributed to two types of nonlinear behavior, both 
of which are illustrated in Figure 3. Referring to the 
figure, curve A experiences an inflection point or a maximum. 
Either condition causes a zero or nearly zero slope and a 
convergence failure. This curve is typical of softening 
nonlinearity. Curve B, on the other hand, is typical of a 
stiffening nonlinearity and a convergence failure will occur 
when the slope becomes large. It is necessary, then, to 
examine the load-displacement behavior of the shell to 
determine the cause of convergence failure. A detailed 
discussion of the theory and implementation of SATANS can 
be found on pages 8 through 29 and in Appendices B, C and D 
of Ref. 3, and on pages 8 through 28, pages 56 through 61 
and page 100 of Ref. 4. 

C. USER PREPARED SUBROUTINES 

Five small subroutines describing the shell geometry, 
the stiffnesses, the loading (physical and/or thermal) , 
and the initial conditions (in the dynamic case) must be 
prepared by the user. Complete descriptions for the 
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LOAD 




DISPLACEMENT 



NONLINEAR BEHAVIOR AFFECTING CONVERGENCE 

FIGURE 3 
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preparation of these subroutines can be found on pages 34 
through 38 of Ref. 4. 
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III. CONSTRUCTION OF SATANS-I 



Stein /Ref. 107" observes that a basic impediment to the 
widespread use of existing computer programs is the difficult 
and laborious process of assimilation of a program's basic 
theory, its implementation and its use by a prospective 
researcher. He also warns the 'unwary analyst' of the dan- 
gers inherent in computer output, which is a 'blizzard of 
numbers.' As a consequence, SATANS is continually being 
modified in an attempt to make it more usable. The most 
recent modifications made for this purpose are described in 
this section. This new version of SATANS is called SATANS-I. 

A. PACKAGING 

SATANS has been reorganized to allow 99% of the code to 
be prepackaged for repeated use by researchers. Compiled 
versions of all subroutines and main controlling routines 
can now be placed on a disk, in a data cell or on a magnetic 
tape, allowing the user to 'submit' only the user written 
subroutines and input data. A complete listing of this 
storable package is given in Appendix E. 

B. PLOTTING ROUTINE 

A great advantage is gained if the computer output can 
be presented in graphical format. Such plots, however, are 
often generated using mechanical or photo-electronic devices 
which are not universally usable. Consequently, SATANS has 
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been enhanced with an equipment-independent plot package 
utilizing self-contained subroutines and the high-speed 
line printer. 

The main program listed in Ref. 4 could not be modified 
to incorporate the logic and control for a graphical ability 
as it then exceeded the capability of the FORTRAN-H compil- 
er. Therefore, the main program of Ref. 4 has been replaced 
by a new main program and two subroutines described as fol- 
lows : 

SATANS - This is the new 'main. ' It reads all data 
cards, prepares the title page of output and passes control 
to either of two basic controlling subroutines, 

STATIC - This is the controlling subroutine for a static 
analysis. It conducts initialization, establishes geometry, 
stiffness, loading and boundary conditions, controls plot- 
ting of input information and controls the solution proce- 
dure. 

DYNAMC - This is the controlling subroutine for a dynam- 
ic analysis and it operates in a fashion similar to subrou- 
tine STATIC. 

The graphical ability is obtained by the introduction of 
two controlling subroutines, PL0T1 and PL0T2, four plotting 
subroutines, PLOTIT, SCALIT, ROUND and DRAWIT,^" and modifi- 
cation of subroutine OUTPUT for the selection of output data 
for plotting. A single data card with request information 

■'"The four plotting subroutines are slightly modified 
versions of those contained in Ref. 11 under the listing 
"JS UTPLOT . " 
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is required to generate the plots. The plotting routines 
will handle a maximum of 100 stations. This routine is 
included in the listing given in Appendix E. 

C. MISCELLANEOUS ALTERATIONS 

The iteration lower limit has been raised to two itera- 
tions. This change was introduced because it was noted 
that immediately after a load step size reduction, several 
consecutive solutions were obtained in one iteration, a 
completely linear process. This was then followed by non- 
convergence and, in some instances, premature program termi- 
nation. Setting a minimum of two iterations automatically 
introduces nonlinearity into each load step and appears to 
solve this problem, as illustrated in Figure 4. 

A quadratic extrapolation technique was investigated 
using a static analysis of the highly nonlinear shallow 
spherical cap. No gain in solution accuracy was seen and 
a price was paid in program execution time due to approxi- 
mately 10% greater total iterations required, not to mention 
the increased storage space requirement. This phenomenon is 
exactly as observed by Stricklin et al /Ref. 12 / , hence the 
linear extrapolation technique has been retained. 

Since boundary loads are introduced through the boundary 
conditions and since the program considers the boundary 
conditions to be the same in every harmonic, changes indi- 
cated in Appendix B of Ref. 5 are permanently incorporated 
into the program. This alteration causes the applied 
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Load 










Deflection 

PREMATURE PROGRAM TERMINATION 
FIGURE 4 
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loading at a boundary to be taken in the zeroeth harmonic 
only, i.e., when n=0. 

All 'COMMON' statements have been ordered and separated 
and some slightly modified. This alteration allows easier 
reading and correlation of FORTRAN variables with specific 
subroutines . 

The output ability has been expanded to allow printed 
data (and plots) at up to 3 6 vice six meridians. 

In the case of a shell with an initial (or a final) pole 
determination of the shear force at the pole in the n=l har- 
monic has been deemed an inconsequential and lengthy computa 
tion and has been removed. Appropriate changes have been 
made in subroutines INLPOL, FNLPOL and OUTPUT. 

D. USE OF SATANS-I 

Instructions for preparation of input da^ta, cards for 
SATANS-I are contained in Appendix C and the initiating main 
program and user written subroutines appear in Appendix D. 
Appendices D and E contain a complete listing of the 
computer program, SATANS-I. 

SATANS-I requires 206 k-bytes of main core when run on 
an IBM-360/67 computer under OS/MVT and compiled by the 
FORTRAN-H compiler at optimization level 2. Execution times 
are the same as for SATANS unless plotting is requested, 
and then is increased approximately 0.01 seconds for each 
100 point plot generated. 
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IV. CONSTRUCTION OF SATANS-II 



As discussed previously, a serious impediment to the 
consideration of initial imperfections by Dantone /Ref. 5/ 
was the requirement that the imperfections be symmetric 
about a datum meridian. This restriction results from the 
fact that all dependent variables were expanded in a single 
sine or cosine series. The analysis presented in this 
section removes this restriction by introducing the use of 
full trigonometric series expansions for the dependent 
variables, and then modifies the resulting new relations 
where required for the introduction of initial imperfections. 
The computer program that incorporates these changes is a 
modification to SATANS-I, called SATANS-II. 

A. INTRODUCTION OF FULL TRIGONOMETRIC SERIES EXPANSIONS 
The original assumption by Ball /Ref. 3 ~J was that a 
normal pressure load could be expressed in a Fourier cosine 
series as 

oo 

q = (o o h Q /E o ) l p (n) Cos n6 (1) 

o 

where a Q = the reference stress level, 
h Q = the reference thickness, 

E q = the reference Young's modulus, 
n = the Fourier index, or harmonic, 
and p = the nondimensional Fourier coefficient of the 
normal pressure load. The normal displacement W 
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corresponding to this pressure was taken as 



• uu 

W = (aa 0 /E 0 ) l w (n) Cos n6 (2) 

o 



where a is the shell characteristic length and w is the 
nondimensional series coefficient of the normal displace- 
ment. For a completely arbitrary case, representing q and 
W with the complete Fourier series gives 

q = (°o h o /E o^ ( l P l(n) Cos n6 + l p 2 (n) Sin n0) (3) 

o 1 



and 



W = (ao o /E Q ) ( £ w 1 ^ n ^Cos n0 + £ w 2 ^ n ^Sin n0) 

o 1 



(4) 



where superscript 1 indicates a previously existing coef- 
ficient and superscript 2 indicates a new coefficient. 
Similar series expansions can be written for the other 
pressures, displacements, etc. 

The nonlinear terms associated with the full series are 
illustrated by 4> 2 , 



4> 2 
s 



(°<> /E o) 



f 

2 



( l * ,(nl Cos no} 2 + ( l ^ 2 < n) sin n6) 2 
^ s s 



+ 2( £ <(>' ( n )cos n0} ( £ if 2 ^ n ^Sin n0) ) 



(5) 



where the <{) are the nondimensional series coefficients of 
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the meridional rotation, given as 



<j) 1 w 1 + u) u 1 (n) 

s s 

^< n > = (-„=< n >)' + Ms u ! <n) 



( 6 ) 



Originally, in Refs. 3 and 4, the nonlinear term 4>| was 
represented as 

<t> 2 = (o 0 /E 0 ) 2 ( l ‘J* g (n) Cos n0) 2 = a 0 /E 0 l B* (n) Cos n0 (7) 

0 o 

since 

Cos (n0 ) Cos (m0 ) = %Cos(n-m)0 + *5Cos(n+m)0 

/ n \ 

The 0 1 ' ' term contains the products of all the <p 1 s whose 
s s 

harmonic sum or difference equals n. There are now two 
additional terms in equation (5). These can be represented 
as 

OO QQ 

(°o /E o) 2 ( l (n) Sin n0) 2 = a Q /E 0 [ 6 2(n) Cos n0 (8) 
i i 

and 

OO oo 

2(° 0 /E 0 ) J ( l *‘ (n) Cos ne)( l Sin n6) 

0 1 

(9) 

= o q /E 0 l 0g 2 (n) Sin n0 
1 

where superscript 1 2 indicates the product of a previously 
existing coefficient and a new coefficient. 
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Thus , 



w w 

= a o/ E o ( I Bg^ n ^Cos n6 + l Bg^ n ^Cos n0 



+ l B* 2(n) Sin n0) 



( 10 ) 



Similarly, the nonlinear term $ 2 (and also $ 2 since it is of 
the same form) is given by 



*§ - K/EoW ( l *i <n, Sin ne) ! + ( J |J<"'co S n0) 



l 1 



+ 2( I *‘ (n) Sin ne}( l ^ Cn) C os n 0 ) ' 



( 11 ) 



which can be represented as 



$ 2 = a 0 /E 0 ( l B* (n) Cos ne + l B* (n) Cos n0 



+ l Bq 2 (n) Sin n0 ) 



( 12 ) 



The nonlinear term $ $ is given by 

s 0 



V 0 = K /E o) 2 ( l 4>g (n) Cos n0j ( l (f g Sin n0) 



+ ( l 4>* (n) Cos n0) ( l <f> 2(n) C os n0) 



+ ( l 4>s (n) Sin n0) ( £ <£* (n) Sin n0) 



( 13 ) 



+ ( l 4>^, (n) Sin n0) ( £ cf> 2 (n) Cos n0) 
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which can be represented as 



OO 



00 




l 2 J n) Cos ne 
se 



(14) 



Following this scheme, similar expressions can be obtained 
for the remaining nonlinear terms $ s N s , ^qNq, ^s N s0' ^0 N 0 S ' 
<J>N S and $N 0 . 

All of the nonlinear terms fall into one of the three 
categories Cos (k0 ) Cos (Z9 ) , Sin (k0 ) Sin (£0 ) or Cos (k0) Sin (£0) . 
Since these three forms were already accounted for in 
SATANS, there was no need to develop any new routines for 
the full series expansions. It was only necessary to intro- 
duce each new nonlinear term into the appropriate existing 
summation in subroutines PHIBET and TEAETA. 

B. INTRODUCTION OF INITIAL IMPERFECTIONS 

The initial imperfections only influence the nonlinear 
terms in the equations. The strain-displacement relation- 
ships used in this analysis are 



e = U' + W/R + h[$ 2 + $ 2 ) 
s s v s 1 



z 



0 



= v‘/r + r"U/r + W/R Q + ^(S 2 + 4> 2 ) 



( 15 ) 



z 



S0 



= h[v' + U*/r - r'V/r + 4> s 4> Q ) 
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where 



$ = -vT + U/R 

s s 

$ 0 = -W/r + V/R 0 (16) 

<i> = Js(v" + r"V/r - U*/r) 



and where e , e A , e represent the reference surface 
s' 0' s9 

strains, $ , 4> , $ represent the reference surface rotations, 

S 0 

and superscript primes and dots represent partial differen- 
tiation with respect to s and 0, respectively. 

Consider a shell whose geometry differs slightly from a 
shell of revolution. For this imperfect shell, a deviation 
from a specific shape occurs in the form of a transverse 
displacement, W. This W produces no strain. Defining W* 
as a total displacement, 

W* = W + W (17) 



and replacing W in equations (15) with W* of equation (17) 
gives 



e + c = 
s s 



e e + £ e = 



e se + e s 6 = 



u' + w*/r + f < f | * 2 + $ 2 ) 
s s 

V*/r + r'U/r + V7*/R Q + + $ 2 ) 

h(v' + u'/r - r'V/r + $*$*) 
v s o 



where 



= -w*" + U/R 
s ' s 

= -w* */r + V/R 0 



(18) 



( 19 ) 
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and e , e A and e A are the pseudo strains due to W only, 
s u s u 

When U = V = W = 0, then W* = W, and £ s = £ 0 = e s0 = 0. 
Substituting these relationships into equation (18) , solving 
for Eg, Tq, ~i S Q and placing these terms in equations (16) 
lead to 

e = U' + W*/R + h 0* 2 + ^ 2 ) - W/R - h(W') 2 
s s ' s s (20) 

e 0 = V* /r + r'U/r + W*/R 0 + + O " W/R 0 " 3 2 (W*) 2 /r 

e s0 = + U*/r - r'V/r + 4>*4>* - W'W*/r) 

Expansion of equations (20) making use of equations (16) , 

(17) and (19) lead to 

£ = U' + W/R + kU 2 + <S> 2 ) - {W"$ } 

s s s s 

e 0 = V*/r + r'U/r + W/R 0 + h {<$> 2 + 4> 2 ) - {W*<I> 0 /r} (21) 

e s0 = h(v' + U*/r - r'V/r + 4> s <I> 0 - {W"<D 0 + W'^/r}) 



where the terms enclosed in braces are the newly generated 
terms due to the imperfections, and 



$* = $ - {w"} 

s s 



$* = $9 _ { W */r} 



( 22 ) 



The presence of imperfections has no effect on the changes 
of curvatures since the curvature-displacement relations are 
linear . 
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Expressing the imperfection partial derivatives -W' and 
-W*/r as imperfection rotations $ s and $ 0 , the strain- 
displacement equations (21) can be expressed as 



where the terms in braces are the newly generated terms. 
Imperfection rotations also affect the nonlinear term, $ S N S 
(and the others listed after equation (14) ). For these 

terms, the reference surface rotations $ s and <i> 0 must be 
replaced by 3>* and 3>* of equations (22) . 

S Q 

To illustrate the procedure which is used in the follow- 
ing section to introduce the initial imperfection rotations, 
consider the single series expansion, equation (7), repeated 
here , 



Modifying this equation to include the effect of initial 
imperfections indicated by equations (23) results in 



e„ = U' + W/R + hi + {2$ $ } + $ 2 ) 

& b v s s s J 

e Q = V*/r + r"U/r + W/R Q + h ($ 0 + + $ 2 } (23) 

e s6 = H{V' + U*/r - $ s $ 0 + {$ 0 ? s + 4 > s 4T q }) 




’ (n) Cos n6 (7) 




Cos ne) 2 



oo 



CO 



+ 2 ( l ip 1 ^ Cos n0) ( 
o s 




( 24 ) 
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A redefinition of 6^ to include the cross-product term is 
now required, along with redefinitions of the other 6 terms, 
and the n terms in the other nonlinear term expressions. A 
redefinition of these terms alone is sufficient for the 
introduction of initial imperfections into a SATANS analysis 
where only the single series expansions are used. This was 
the thrust of the thesis by Dantone. 

C. COMBINATION OF FULL TRIGONOMETRIC SERIES EXPANSIONS AND 
INITIAL IMPERFECTIONS 

Consider that for every meridional rotation $ , where 

*S = (°o/ e o) 2 ( l <t>s Cn) Cos net? <J> 2(n) Sin ne) (.25) 

o 1 

there exists an initial meridional imperfection rotation 
4> s , defined as 

oo oo 

¥ = (o 0 /E 0 ) 2 ( l ^ (n) Cos'ne + l 4>" 2 (n) Sin n0) (26) 

The nonlinear term previously given by equation (5) , now 

becomes (4> 2 + 2$ s $ s ) , according to equation (23) , where $ 2 
is given by equation (5) and is the product of equations 

(25) and (26) , as illustrated for the single series expan- 
sion in equation (24). All the nonlinear terms of equations 
(15) have been converted to those in equation (23) for the 
full series expansions and the resulting relations are 
given in Appendix A. The definitions of the 3 terms (given 
previously in Ref. 3) have been modified to include the 
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effects of both full series expansions and initial imperfec- 
tion presence. 

The remaining nonlinear terms (those listed after equa- 
tion (14) ) also require modification but are of simpler 

form since only the first power of the rotation is involved. 
A redefinition of the n terms to include the effects of full 
series expansions and appropriate initial imperfection rota- 
tions has also been made. 

In this manner, the validity of the strain-displacement 
equations (equations (5) of Ref. 3 and equations (15) of 
this report) and the equilibrium equations (equations (8) 
of Ref. 3) is retained. Introduction of full series 
expansions and initial imperfections in the nonlinear terms 
was accomplished by a redefinition of the 8 and n terms, 
which are the nondimensional coefficients for the nonlinear 
terms in the strain-displacement and equilibrium equations, 
respectively . 

D. SUBROUTINE MODIFICATIONS 

Introduction of the full series into the program was 
implemented by use of negative index numbers for the new 
sine or cosine portions of the expansions (the terms with 
the superscript 2 ) . With the requirement that the axisym- 
metric harmonic (n=0) be included, the program was expanded 
to handle a maximum of 25 modes, i.e., the zeroeth harmonic 
(n=0) plus twelve other harmonics (positive and negative) . 
Each positive index and each negative index will be referred 
to as a mode, i.e., the harmonic n=5 has two modes, -5 



37 



and +5. A typical correspondence between harmonic, mode 
and Fourier index is given in Table I. Harmonic zero (n=0) 
and any combination of up to twelve harmonics may be 
included in the analysis. The positive values of NN in 
Table I are the modes analyzed in SATANS and SATANS-I. The 
negative values of NN are those modes added to incorporate 
the full series expansions in SATANS-II. Note that n=0 is 
the first mode, and that the negative value is always the 
first mode of a pair of modes. This is a requirement. 



TABLE I 

HARMONIC, MODE AND FOURIER INDEX CORRESPONDENCE 



HARMONIC 


MODE 


FOURIER 


FORTRAN ASSIGNMENT 


n 




INDEX 


STATEMENT IN SUBROUTINE 
PLOAD OR TLOAD 


0 


1 


0 


NN ( 1 ) = 0 




2 


-1 


NN (2) =-l 


1 


3 


+1 


NN ( 3 ) =+l 


o 


4 


-2 


NN (4) =-2 


£ 


5 


+ 2 


NN ( 5) =+2 


r 


6 


-5 


NN ( 6) = -5 


D 


7 


+ 5 


NN ( 7 ) =+5 


• 

• 

• 

18 


• 

• 

• 

24 


• 

• 

• 

-18 


• 

• 

• 

NN (24) =-18 


25 


+18 


NN (25) =+18 
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Seven subroutines of SATANS-I required changes to 
incorporate the full trigonometric series expansions and 
the initial imperfections. 

Subroutine PHIBET calculates the <f> ' s and carries out 
the summation procedures for calculating the 3 nonlinear 
terms. Since extensive modification was required, the 
subroutine was altered such that those sections handling 
the full trigonometric expansions and imperfections are 
entered only if required, i.e., if a full series analysis 
is necessary. 

Subroutine TEAETA calculates the inplane forces and 
the n nonlinear terms. It was also extensively modified, 
again, utilizing new sections only if required. 

The subroutine MODES, a highly complex procedure for 
determining modal coupling required minimal alteration 
since each unique mode number is stored at every, or every 
other storage location, depending on use of SATANS-I or 
SATANS-II, respectively. Each unique harmonic is accessed, 
then, by an appropriate step size change in the controlling 
loops . 

The subroutine OUTPUT required the addition of a section 
to include the new half of the trigonometric series in the 
generation of a summed solution at a specified meridian. 

The subroutines POLE, INLPOL and FNLPOL compute the 
solution at an initial or a final pole. Due to the compu- 
tation of pole conditions based on n, which now has positive 
and negative values, these routines were modified 
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accordingly. Effectively, the only change in the equations 
presented in Appendix D of Ref. 3 involve a replacement of 
n by -n where appropriate. 

A complete listing of these seven modified subroutines 
is given in Appendix F. 

E. SATANS-II 

To modify SATANS-I to handle completely arbitrary 
loading and initial imperfections, the following steps 
must be performed. 

1. Replace subroutines PHIBET, TEAETA, MODES, OUTPUT, 
POLE, INLPOL and FNLPOL with those in Appendix F. 

2. Replace the affected COMMON statements in Appendices 
D and E with those listed in Appendix F. 

3. Add subroutine IMPERF, also listed in Appendix F. 
This subroutine contains the values of the imperfection 
rotation coefficients, <j) and <J>g. Instructions for actual 
preparation of these coefficients will be discussed in de- 
tail in section V. If imperfections are not to be consider- 
ed, the <jT g and are never accessed. 

4. In Appendix E, make the following additions after 
the card indicated: 

Subroutine After Card Number Insert 



DYNAMC 



00002270 



CALL IMPERF 



STATIC 



00002230 



CALL IMPERF 



PMATRX 



00001100 IF (N (MN) .LT . 0) P (2 , 1 , IJ) =-l . 0 



PMATRX 



00002000 IF(N(MN) .LT. 0) ZFP1 ( 2 , 2 ) =1 . 0 



SATANS 



00000380 READ (5,108) JUMP, MPERFS 
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SATANS 



00001400 



108 FORMAT (21 5) 



5. In Appendix E, modify the following cards: 
Subroutine Card Number Should Read 

DYNAMC 00000800 DIMENSION VBAR ( 25) ,AVB ( 25) 

PMATRX 00000900 NN=IABS (N (MN) ) 

PMATRX 00001640 NN=IABS (N (MN) ) 

6. In Appendix C, modify the data preparation instruc- 
tions by adding a data card after card #4 according to the 
following format: 



Card 


Column 


Format 


Item 


Interpretation 


4.1 


1-5 


15 


JUMP 


1 = single series expansion 

2 = double series expansion 


4.1 


6-10 


15 


MPERFS 


0 = no imperfections 



1 = imperfections 



An initial imperfection solution requires the use of full 
series in the analysis. Specification of the parameter 
JUMP indicates whether the solution will be as in SATANS-I 
or will utilize the full trigonometric series expansions of 
SATANS-II. Specification of the parameter MPERFS indicates 
whether the analysis will include initial imperfections. 

Note that if MPERFS equals 1, JUMP must equal 2, but if 
JUMP equals 2, MPERFS may be either 1 or 0 . 

7. Add imperfection data just prior to the last card of 
the data deck. This data is to be placed in arrays PHIXB 
and PHITB. 
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F. USE OF SATANS-II 

After completing the changes in section IV. E, SATANS-II 
exists as a program capable of performing all the functions 
of SATANS-I with the following restrictions: 

1. The number of meridional stations cannot exceed 
100, (KMAX < 101) . 

2. The number of modes cannot exceed 25, (MNMAX < 26) . 

3. The product of meridional stations and modes cannot 
exceed 951, (KMAX* MNMAX < 952) . 

4. Plotting ability is not available for the imperfec- 
tion input data. 

Specification of harmonics to be considered (accomplish- 
ed in subroutine PLOAD or TLOAD) must include the axisym- 
metric harmonic (n=0) and both negative and positive entries 
for any other harmonic, as indicated in Table I. 

The instructions for the preparation of the input data 
are given in Appendix C. SATANS-II requires 426 k-bytes of 
main core on an IBM-360/67 computer when run under OS/MVT 
and when compiled by the FORTRAN-H compiler at optimization 
level 2. Execution times vary widely from four minutes for 
a 35 station and 3 mode analysis, up to 30 - 60 minutes for 
a 35 station and 25 mode analysis. 

G. VALIDITY OF SATANS-II 

SATANS-II was checked for its ability to obtain a cor- 
rect solution using the full series expansions. This check 
was made using a clamped cylindrical shell subjected to 
pressure loadings as described in Figure 5. A solution was 



42 






4- = SOLUTION COMPARISON POINTS 

SATANS-IE CHECK LOADINGS 
FIGURE 5 
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obtained for the loading shown in Figure 5A with SATANS-I 
and the summed solution was printed at meridian 0=0°. The 
cylinder was then analyzed by SATANS-II for the following 
three orientations: 

1. Loading defined as in Figure 5A - a symmetric load. 

2. Loading defined as in Figure 5B - an anti-symmetric 
load. In this case, the summed solution was printed at 
meridian 0 = 90°. 

3. Loading defined as in Figure 5C - a general load. 

In this case, the summed solution was printed at meridian 
0 = 45°. 

Load definitions for the four check runs appear in 
Table II. In each case, the solution obtained by SATANS-II 



TABLE II 

FOURIER COEFFICIENTS FOR CHECK LOADINGS 



MODE 


FOURIER COEFFICIENTS OF APPLIED LOADING 


SYMMETRIC 


ANTI-SYMMETRIC 


NONSYMMETRIC 


0 


0.0 


0.0 


0.0 


-1 


0 . 0 


-63.66 


-45.015 


+1 


-63.66 


0.0 


-45.015 


-3 


0.0 


+21.22 


+15.005 


+3 


+21.22 


0.0 


-15.005 


-5 


0.0 


-12.73 


+9.003 


+ 5 


-12.73 


0.0 


+9.003 


-7 


0.0 


+9.095 


-6.431 


+7 


+9.095 


0.0 


+6.431 



Note: the coefficients are based on p Q = 50 psi. 
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2 

agreed with the SATANS-I solution through four significant 
digits. 



This reference to solution includes the load-deflection 
history, the ultimate nonconvergence level, and the displace- 
ment, force and moment distributions. 



45 



V. IMPERFECT SHELL ANALYSIS 



A. PREPARATION OF THE IMPERFECTION DATA 

With the developments of section IV implemented, 

SATANS-II requires only the Fourier coefficients of the 

imperfection data to conduct a nonlinear analysis. Recall 

from equation (26) that these Fourier coefficients, <f> s and 

are determined by knowledge of the actual imperfection 

rotations, ¥ and ¥ fl . Recall also from equation (22) that 
S o 

the imperfection rotations are given as 



$> = -vr 

s 

¥ 0 = -W*/r 



(27) 



This means that the imperfection rotation Fourier coeffi- 
cients can be determined by knowledge of the initial normal 
deflection deviation, W. This imperfection information can 
be obtained from a mathematical description, a graphical 
plot or from actual experimental measurement. A simplified 
data generation flow chart is presented in Figure 6. 

1 . Data Generation by Mathematical Description 

In a note on imperfect shallow spherical cap buck- 
ling, Kao /Ref. 137 analyzes a cap subjected to a mathemat- 
ically defined dimple imperfection existing over one-quarter 
of the shell. The dimple is given by 
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SIMPLIFIED IMPERFECTION DATA GENERATION 

SCHEME 

FIGURE 6 
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W = 64 (s/l) 3 (1-s/l) 3 h6Cos 20, -45°<6< 45° 



( 28 ) 



W = 0, 45 O <0<325° 

where L is the total length of the shell meridian and where 
s/L ranges from 0 at the pole to 1 at the edge, and 6 is an 
intensity factor, 0 <= 6 <= 1. The geometry of this problem 
is shown in Figure 7. Kao used a two-dimensional finite 
difference mesh of eight meridional and 16 circumferential 
stations in his analysis. A numerical model of, for example, 
35 meridional stations and 25 modes would yield a SATANS-II 
solution of considerable accuracy. The discrete initial 
normal deflection deviations W at each mesh point need not 
be computed, as partial differentiations of equation (28) 
give expressions directly solvable for discrete values of 
4> and $ for each s and 9. Imperfection rotations at loca- 

S u 

tions exterior to the affected area (-45° < 9 < 45°) are 

zero. To illustrate, the meridional imperfection rotations 

that exist along the line ABCDEA of Figure 7 would appear 

as curve A of Figure 8. The discrete values generating 

this curve can easily be expanded into a full Fourier series 

( n ) 

and the resulting coefficients d> ^ inserted into SATANS-II 

s 

for the analysis. 

The imperfection of Figure 7 is symmetric about the 
meridian 0=0°. As an example of a more complex imperfec- 
tion, another dimple of the same form could be located in 
the region defined by 100° < 0 < 190°. In this asymmetric 
case, the meridional imperfection rotations along line 
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DIMPLE IMPERFECTION GIVEN BY 
W/h = 64X 3 (|-X? cf Cos 20 
WHERE X = yL 

MATHEMATICAL IMPERFECTION GENERATION 

FIGURE 7 
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MERIDIONAL IMPERFECTION ROTATION 

REPRESENTATION 

FIGURE 8 
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ABCDEA of Figure 7 would now appear as curve B of Figure 8. 
The discrete values generating this curve can also be easily 
expanded into the full trigonometric series and its coeffi- 
cients inserted into SATANS-II. 

2 . Data Generation by Graphical Plots 

Graphical representations of W can be easily con- 
structed and the patterns are near limitless in design or 
location. Theoretical or actual creases, furrows or dents 
of any shape need only be graphed, the discrete imperfection 
displacements determined, the discrete imperfection rota- 
tions ¥ and $ computed via a finite differencing tech- 
s 0 

nique, the Fourier expansions performed and the resulting 
coefficients inserted into SATANS-II. 

3 . Data Generation by Experimental Measurement 
Arbocz and Babcock /Ref. 67 have developed a means 

of mapping the topography of circular cylindrical shells, 
and can thus describe the deviations of an actual shell. 

The procedure, described in detail in Ref. 6, produces 
imperfection data in the form of nondimensional normal 
displacements from a perfectly circular cylinder. For the 
sample shells considered here, 49 data points were taken 
around each circumference and 29 equally-spaced circumfer- 
ential traces were made. The first and last trace were 
0.64 inches from the ends of the shell and the spacing 
between traces was 0.24 inches. Pertinent geometrical and 
material data concerning the shells of interest here is 
contained in Table III. 
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TABLE III 



IMPERFECT SHELL OVERALL CHARACTERISTICS 



SHELL DESIGNATION 


A- 13 


A-14 


Radius 


4.00" 


4.00" 


Length 


8.00" 


8.00" 


Young's Modulus x 10 6 , E 


15.10 psi 


15.80 psi 


Wall thickness, h 


0.00444" 


0.00437" 


Number of axial traces 


29 


29 


Axial increment, A 


CNJ 

• 

o 

1 1 


0.24" 


Number of data points 


1421 


1421 



To prepare the imperfection data on shell A-14 for 
use by SATANS-II, the nondimensional deviations were first 
dimensionalized by the wall thickness. To retain the 
original axial increment of 0.24 inches between circumfer- 
ential traces and to make the data compatible with a 35- 
station meridian, the shell length was increased by 0.08 
inches at each end, or a total of 2% (See Figure 9). This 
required three additional traces at each end of the shell 
corresponding to stations 1, 2 and 3, and 33, 34 and 35. 

Data for the additional six circumferential traces were 
generated via the modified linear extrapolation procedure 
illustrated in Figure 10. At each of the 1,715 data points, 
the meridional and circumferential rotations were computed 
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8. 00 inches 




CIRCUMFERENTIAL 
TRACE NUMBERS 
(REF. 6) 



MERIDIONAL 
STATION NUMBERS 
(SATANS-H) 

NOTE: Extrapolation required for dato 
at primed stations 



CIRCUMFERENTIAL TRACE NUMBER AND 
MERIDIONAL STATION NUMBER CORRESPONDENCE 

FIGURE 9 
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I 

I 



• = KNOWN DATA POINT (REF. 6) 

B = TW 0 STATION AVERAGE 
& - EXTRAPOLATED DATA POINT 

Wjj= IMPERFECTION DISPLACEMENT 
J OF THE ith STATION AND THE 
j th MERIDIAN 

NOTE: NUMBERS CORRESPOND TO 
SATANS -H. MERIDIONAL STATIONS 



DATA POINT EXTRAPOLATION TECHNIQUE FOR W 

FIGURE 10 
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using central finite differencing in the appropriate direc- 
tion. After nondimensionalization according to the 
SATANS-II scheme, each circumferential vector of 48 
meridional and 48 circumferential imperfection rotations 
was expanded via a periodogram into a Fourier series accord- 
ing to 

C E o/ a oH s = + I 4 ^ (n) Cos n9 + ^ (n) Sin n 6 (29) 

and 

( E o / a o ^0 = ^ 0 °^ + I ^^^Cos n0 + ^ 4>g^ n ^Sin n0 (30) 

At this point, the Fourier coefficients <p and 4> are of 

s 0 

proper form for direct insertion into SATANS-II. The 
computer program which performed the actual imperfection 
data processing for shell A-14 is given in Appendix G. 

4 . Selection of Important Harmonics 

Since the number of meridional stations was 35, the 

3 

program maximum of 25 modes can be used in the analysis. 

Note, however, that in equations (29) and (30) , there are 
49 modes describing the imperfections. Therefore, a selec- 
tion of a maximum of 25 imperfection modes must be made. 

The harmonics to be considered in the analysis should in- 
clude those critical to the shell and those significant 
in the imperfection. The description of the harmonic 

^Recall the restrictions given m section IV. F. 
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selection process for shell A-14 is given in a later 
section. 



B. BOUNDARY CONDITIONS 

The testing configuration used by Arbocz and Babcock 
on the axially loaded cylinder is shown in Figure 11. The 
ends of the shell are attached to an end ring and a load 
cell with cerrelow, a low melting-point alloy. The ends of 
the shell are believed to be clamped so that 

W = $ s = V = 0 (31) 

The axial load on the test specimen was applied by means of 
four screws located at the corners of the square end plates. 
At each load increment, the screws were adjusted to give as 
uniform a circumferential load distribution as possible. 

In all versions of SATANS, the boundary conditions are 
represented by the matrix equation 



f 

N 

s 




' 'i 

U 


N 

s 


+ / A _7 


V 


Q s 




w 


. *s . 




M 

sj 



U) 



(32) 



Boundary conditions can be specified on N s or U, N s g or V, 

Q s or W, and £> s or M s by prescribing input values for the 
matrices / ft_/ and / — A J . There is one equation (32) for 
each boundary of the shell. In their present form, SATANS-I 
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CYLINDRICAL SHELL TESTING CONFIGURATION 

FIGURE II 
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and SATANS-II assume the column matrix {£} as zero for all 
harmonics except the axisymmetric harmonic. 

As was mentioned earlier, boundary condition discrepancy 
between mathematical model and physical experiment is cited 
as a major contributor to thin shell buckling load disagree- 
ment. For this reason, three interpretations of Figure 11 
were made. The cylinder was represented as experiencing an 
axisymmetrically applied force N g on both ends (condition A) , 
an axisymmetrically applied force N g on one end and a 
meridional deflection restraint on U at the other (condi- 
tion B) , and an axisymmetric incremental end shortening U 
at one end with meridional deflection U restrained at the 
other (condition C).^ These three conditions were investi- 
gated on the clamped cylinder, and on a simply supported 
cylinder for further comparisons and validity checks. 
Mathematical assignments for boundary conditions A, B and 
C are shown in Table IV. 

C. ANALYSIS 

1 . Buckling of the Thin Perf ectly-Circular Cylinder 
A necessary part of this attempt to develop a pro- 
gram for predicting the buckling load of imperfect shells is 
the establishment of the program's ability to predict the 
buckling load of perfect shells. Consequently, this section 
compares the results of a SATANS-II analysis of the thin 



It is believed that boundary condition C best repre- 
sents the actual boundary conditions in Figure 11. 
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BOUNDARY CONDITION SPECIFICATIONS 
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Note: When or U v ' are specified nonzero, they are set to zero for 



V 



si 



cylinder without imperfections with the bifurcation buckling 
analysis presented by Almroth /Ref. 1£/. Almroth's study 
deals with the influence of edge constraints on the critical 
load of cylindrical shells of the type used by Arbocz and 
Babcock. Equation (33) is given by Almroth for the axisym- 
metric displacement shape of the simply supported cylinder 
(boundary condition A) prior to buckling. 



W 



r (2vyN g + p) (l + A X S in (a ! s) Sinh (a j s) 



(33) 



+ A 2 Cos (a 2 s) Cosh (a 2 s) ) 



where 



a l - k( 1 + Ng/ 

a 2 - k( 1 - Nj* 5 
k =(2y)-*S 

Y = h/r/(l2(l - v 2 ))** 



N 



s 



= -N s /(2yEh) 



p = (pr) / (Eh) 
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Ai 



A 2 



({1 - Ng J^Sin (a j ) Sinh (a 2 ) + N g Cos (a t ) Cosh (a 2 ) ) 
(l - N 2 ) ^(Cosh 2 (a 2 ) + Sin 2 (a 2 ) ) 

({1 - N|} 2 Cos(ai)CoshCa 2 ) - NgSin (a x ) Sinh (a 2 ) ) 
(1 - N 2 ) ^(cosh 2 (cx 2 ) - Sin 2 (a 2 ) } 



aj = (a 1 L)/(2r) 
a 2 = (a 2 L) / (2r) 



and where N g is the axial load, p is a uniform normal pres- 
sure, h is the thickness, E is Young's Modulus, v is 
Poisson's ratio, L is the length and r is the radius. This 
solution is referred to as the consistent prebuckling solu- 
tion, as opposed to the membrane or classical buckling solu- 
tion where W is a constant along the meridian. Equation 
(33) has been solved for shell A-14 and W is plotted as a 
function of s in Figure 12 for an N g = 0.21904. Also shown 
in Figure 12 is a plot of W versus s obtained from SATANS-II 

for the same N for both 35 and 69 meridional station anal- 
s 

yses. Note in Figure 12 that even though only a few sta- 
tions lie in the boundary zone, the trend is towards the 
correct solution. The displacement of the 35 station solu- 
tion, however, is only a rough approximation of the actual 
state . 

Almroth has determined the minimum bifurcation 
buckling load, N Scr , and the corresponding value for n for a 
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MEMBRANE PREBUCKLE DEFLECTION SHAPE, SHELL A-14 



range of parameters L/r and r/h, using both the consistent 
prebuckling displacement W given by equation (33) and the 
membrane or classical prebuckling solution, where W is 
constant, in the Donnell equations. For the simply support- 
ed cylinder with an r/h of the same order of magnitude as 

shell A-14 and an L/r = 0.5, the critical load is N_ =0.863 

s cr 

^scr = ^’^ lb/in) with n=27 for the consistent analysis, and 
N Scr =1.0 (N Scr =46 lb/in) with n=27 for the membrane pre- 
buckling analysis. Thus, the consideration of the axisym- 
metric edge effect, which is analogous to an imperfection, 
reduces the buckling load by 14.1%. Parametric studies 
conducted by Almroth show usage of from 100 to 800 merid- 
ional stations over the half-length of the shell. 

As shown in Ref. 5, use of the membrane prebuckling 
displacement in Sanders' equations for the simply supported 
cylinder simplifies the problem to an eigenvalue problem. 

The buckling mode is 

W = 6 Sin (mus/L) Cos (n0) (34) 

The equations for the matrix elements in the eigenvalue 
problem are given in Appendix B. A plot of the critical 
values of m and n is given in Figure 13. The minimum buck- 
ling load occurs when n=27 and m=15, and corresponds to an 

N s C r of 47 * 5 lb / in * 

SATANS ' ability to predict the minimum bifurcation 
load of a shell for a specified value of n has been demon- 
strated by Stilwell /Ref. 15/- SATANS-II analyses of shell 



63 



40 — 




m 



cr 



CRITICAL B'JCKU.NG MODES FROM EIGENVALUE 
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FIGURE 13 
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A-14 (simple support, boundary condition A) were conducted 
over a wide range of n for both the consistent and membrane 
prebuckle cases. Figures 14 and 15 give plots of the crit- 
ical buckling load versus the number of meridional stations 
for various values of n. The dashed parts of the curves 
are estimated extrapolations. It is clear that as the 
number of stations is increased, SATAN S ' critical buckling 
load decreases and the associated harmonic increases. While 
the n=27 curve predicts the highest buckling load over the 
range of stations used in both cases, extrapolation of 
these curves indicates that the n=27 curve would predict 
the lowest buckling load if the number of stations exceeded 
200 to 300. Note that the membrane prebuckle critical 
buckling curves of Figure 15 seem to flatten much more 
quickly than for the consistent case. This appears to 
confirm the occurence of a reduction in critical buckling 
load due to the presence of an axisymmetric edge restraint, 
as noted by Almroth. Figures 16 and 17, which present plots 
of the critical buckling load versus critical buckling mode 
of shell A-14 for the 35, 69 and 100 station analyses, also 
predict this lowering of buckling load due to edge restraint, 
the 'fanning' of curves in Figure 16 indicate an increase 
in the number of meridional stations will produce a lesser 
buckling load decrease than in Figure 17. These figures 
also predict a lower critical buckling load and a higher 
critical buckling mode with increasing meridional stations. 
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Figures 14 through 17 indicate that, given enough 
stations, SATANS-II can predict a critical buckling load 
and mode for the thin-walled shell that agrees with 
Almroth's prediction. However, for the 35-station analysis, 
which is the analysis to be used for the imperfect shells, 
the critical load of 51.1 lb/in and harmonic of 16 are 
significantly different from the correct values. As a 
result of the coarse mesh used by SATANS-II in the following 
analysis, the critical buckling loads for the imperfect 
shells will be higher than the experimental values. 
Nevertheless, it will be shown that the introduction of a 
totally arbitrary, naturally-occurring imperfection pattern 
on the imperfection-free mathematical model of shell A-14 
results in a significant reduction in critical buckling 
load. This is not only sufficient to demonstrate the 
capability of the program to perform the analyses, but it 
also provides an estimate of the degree of imperfection 
sensitivity of the shell. It will be shown that application 
of this estimate of sensitivity to Almroth's buckling load 
provides an accurate prediction of the buckling load of the 
imperfect shell. 

2 . Introduction of Imperfections 

When prepared as indicated in section V.A.2, the 
imperfection rotation Fourier coefficients are in a form for 
direct introduction into SATANS-II. Figure 18 shows the 
applicability of the periodogram method of expansion of the 
imperfection data according to equations (29) and (30). 
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Curve A is the graph of a solution of equation (30) pre- 
senting the circumferential imperfection rotation at 
station 10 plotted against 0. It is seen to pass smoothly 
through all the data points. 

Only the zeroeth harmonic (n=0) and twelve others 
can be accommodated by SATANS-II. In the case where the 
shell response and the imperfection are limited to a few 
major harmonics, all important harmonics can be included in 
the analysis. For example, the X=6 shallow spherical cap 
has only the n=2 harmonic as the major harmonic of response 
and the imperfection curves of Figure 8 require only a few 
Fourier harmonics for adequately accurate representation. 
Also, the few required harmonics for imperfection represen- 
tation are in the vicinity of the n=2 harmonic. This pro- 
blem is will within the limits of SATANS-II. 

The thin-walled cylinder, however, responds in a 
large number of harmonics. This can be seen in Figure 16 
— many harmonics produce critical buckling loads close to 
the value of the minimum critical buckling load. Addition- 
ally, the imperfection requires many harmonics. Thus, as 
mentioned in section V.A.4, a difficult selection of harmon 
ics for use in the analysis must be made. Six harmonics 
were selected on the basis of the shell response and six 
harmonics were selected on the basis of the imperfection 
characteristics. Of course, all 13 selected harmonics will 
contain both the imperfection characteristics and the shell 
response characteristics. 
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For the shell response harmonic selection, the 
critical buckling loads of shell A-14 for specified values 
of n, discussed in the previous section, are presented in 
Table V. Values were assigned to each harmonic according 
to the legend of Table V and six were selected based on 
these assignments. 

For the imperfection harmonic selection, the sums of 
the absolute values of the Fourier coefficients for both 
sine and cosine terms of both meridional and circumferential 
imperfection rotations were tabulated and are presented in 
Table VI. Values were again assigned to each harmonic ac- 
cording to the legend of Table V, and six additional harmon- 
ics were selected. 

The complete harmonic selection process for shell 
A-14 was accomplished by a comparison of the values assigned 
to each harmonic for both the imperfection and the shell. 
Table VII presents this comparison along with the actual 
harmonics selected for the analysis. Note that no overlap 
of important harmonics occurred, as the imperfection re- 
quired harmonics seven and below and the shell required 
harmonics ten and above. The harmonics used with shell A-13 
were n=0 and n=2 , 3 , 4 , 6 , 8 , 9 for the imperfection and n=10,ll, 
12,13,14,15 for the shell response with a simple support, 
and n=13 , 14 , 15, 16 , 17 , 18 for the shell response with a fixed 
support. 

Once the specific harmonics for inclusion in the 
analysis are determined, the appropriate ¥ s and 
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TABLE V 



CRITICAL BUCKLING LOADS OF INDIVIDUAL HARMONICS 



HARMONIC 


FIXED SUPPORT 


SIMPLE SUPPORT 


BUCKLING 
LOAD 
x P Q 


p / p cr 


ASSIGNED 

VALUE 


BUCKLING 

LOAD 

x P o 


p / p cr 


ASSIGNED 

VALUE 


2 


6+ 






6+ 






3 


6+ 






6+ 






4 


5.320 






! 6+ 2 






5 


5.120 






5.720 






6 


4.250 






4 . 080 






7 


2.415 


2.380 




2.445 


2.390 




8 


1.611 


1.550 




1.661 


1.630 




9 


1.221 


1.180 


- 


1.291 


; 1.260 




10 


1.093 


1.050 


* 


1.157 


1.130 




11 


1.088 


1.050 


k 


1.131 


11.110 


- 


12 


1.069 


1.030 


k 


1.085 


1.060 


+ 


13 


1.038 


1.000 


k k 


1.037 


1.014 


k 


14 


1.043 


1.010 


k 


1.029 


1.006 


k 


15 


1.054 


1.020 


* 


1.026 


,1.004 


k 


16 


1.058 


1.020 


•k 


1.022 


1.000 


k k 


17 


1.074 


1.030 


k 


1.030 


1. 008 


k 


18 

19 


1.083 

1.120 


1.040 


★ 


1.037 


1.015 


k 


1.080 


+ 


1.040 


"1.018 


k 


20 


1.200 


1.160 


- 


1.069 


1.050 


-j- 


21 


1.240 


1.190 


- 


1.090 


1.07 0 




22 


1.320 


1.270 




1.130 


1.110 


- 


23 i 


1.340 


1.290 




1.168 


1.144 


- 


24 


1.344 


1.290 




1.200 


1.170 


— 



Notes : 



1 

2 



P Q = 50 psi. 



The trigger in each harmonic is an internal 



pressure of 10 7 psi. 

Each analysis included the indicated harmonic 
plus n=0. 



legend : 



SYMBOL 


VALUE 


k k 


absolutely necessary 


k 


great 


+ 


moderate 


- 


slight 




none 
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TABLE VI 



RELATIVE IMPERFECTION ROTATION COEFFICIENT SIZES , SHELL A-14 



HARMONIC 




c 








4 


>9 




COS 


VALUE 


SIN 


VALUE 


COS 


VALUE 


SIN 


VALUE 


0 


5.32 


* 


0.00 




0.01 




0.00 




1 


3.18 


* 


2.03 


* 


0.34 




0.40 




2 


4.85 


* 


9.08 


* 


9.46 


* 


13.26 


* 


3 


1.50 


+ 


7.34 


* 


13.46 


* 


3.11 


+ 


4 


2.80 


* 


1.74 


+ 


6.28 


+ 


28.54 


* 


5 


1.71 


+ 


0.79 




7.84 


+ 


9.98 


* 


6 


1.33 


+ 


1.48 


+ 


5.20 


+ 


4.07 


+ 


7 


0.56 




1.85 


+ 


11.20 


* 


6.26 


+ 


8 


0. 57 




0.86 


- 


7.88 


+ 


1.92 


- 


9 


0.46 




0.69 




3.75 


- 


2.88 


- 


10 


1.10 


+ 


0.85 


- 


5.60 


+ 


5.99 


- 


11 


0.77 




0.91 


- 


2.17 


- 


1.68 




12 


1.39 


+ 


1.67 


+ 


5.73 


+ 


6.77 


+ 


13 


1.05 


+ 


1.01 


+ 


2.98 


- 


3.38 


+ 


14 


0.77 




0.83 


- 


2.60 


- 


1.55 


+ 


15 


0.91 


- 


0.46 




1.05 


- 


1.26 


- 


16 


1.09 


+ 


0.69 




0.73 




2.17 


- 


17 


0.54 




1.15 


+ 


1.31 


- 


1.17 


- 


18 


0.83 


- 


1.03 


+ 


1.16 


- 


1.46 


- 


19 


0.94 




0.78 




0.48 




1.04 


- 


20 


0.41 




0.81 


- 


0.49 




0.48 




21 


0.42 




0.68 




0.45 




0.19 




22 


0.43 




0.63 




0.31 




0.15 




23 


0.68 




0.58 




0.11 




0.10 




24 


0.32 




0.01 




0.01 




0.01 





Notes: 1. Numbers represent the sums of the absolute values 

of the imperfection rotation coefficients over 
the entire shall. 

2. Value assignments were made according to the 
legend of Table V. 
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TABLE VII 



HARMONIC SELECTION, SHELL A-14 



H 

A 

R 


coefIi- 

CIENT 


<i>e 

COEFFI- 

CIENT 


SHELL 

RESPONSE 


SELECTED 

HARMONICS 


M 




















0 


COS 


SIN 


COS 


SIN 


FIXED 


SIMPLE 


FIXED 


SIMPLE 


N 










SUPPORT 


SUPPORT 


SUPPORT 


SUPPORT 


I 














PRIMARY 


ALTER- 


PRIMARY 


c 
















NATE 




6 


* 












X 


X 


X 


l 


* 


* 










X 


X 


X 


2 


* 


* 


* 


* 






X 


X 


X 


3 


+ 


* 


* 


+ 






X 


X 


X 


4 


* 


+ 


+ 


* 






X 


X 


X 


5 


+ 




+ 


* 






X 


X 


X 


6 


+ 


+ 


+ 


+ 












7 




+ 


* 


+ 






X 


X 


X 


8 

Q 




— 


+ 














10 


+ 


— 


+ 


+ 


* 


— 


X 






11 




- 


- 




* 


- 


X 






12 


+ 


+ 


+ 


+ 


* 


+ 


X 


X 




13 


+ 


+ 


- 


+ 


* * 


* 


X 


X 


X 


14 




- 


- 


- 


* 


* 


X 


X 


X 


15 








- 


* 


* 


X 


X 


X 


16 


+ 






- 


* 


** 




X 


X 


17 




+ 


- 


- 


* 


* 




X 


X 


18 


- 


+ 


- 


- 


* 


* 






X 


19 


- 






- 


+ 


* 


- 






20 




- 






- 


+ 








21 










- 


+ 








22 












- 








23 












- 








24 





















Notes: 1. Harmonic zero is required by SATANS-II. 



2. Value assignments are taken from Tables V and VI. 



3. 'X' indicates a selected harmonic. 
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coefficients (obtained from the periodogram) must be loaded 
into the SATANS-II arrays PHIXB and PHITB (performed in 
subroutine IMPERF) in the following order: 

1. The first KMAX locations contain the <f> 's (or <}> 's) 

s y 

of mode 0. 

2. The second KMAX locations contain the $ s 's (or <Fq's) 
of mode -i. 

3. The third KMAX locations contain the <F 's (or <F's) 

s t) 

of mode +i. 

4. The fourth KMAX locations contain the <p 's (or <fT ' s) 

S u 

of mode - j . 

5. Et cetera, until all modes are entered. 

KMAX is the FORTRAN name for the number of meridional sta- 
tions to be included in the analysis. One special case must 
be considered where the shell has an initial or a final 
pole. In this case, rotations are zero in all but the n=l 
harmonic. If this harmonic is included in the analysis, 
the imperfection rotation coefficients at the first (or 
last) station should all be set to zero except for the 
coefficients in the n=l harmonic. In the case of an initial 
pole, the coefficient associated with the -1 index should be 
placed in the first location of the second set of KMAX loca- 
tions in arrays PHIXB and PHITB, and the coefficient assoc- 
iated with the +1 index should be placed in the f ir st 
location of the third set of KMAX locations. For a final 
pole, the coefficient associated with the -1 index should be 
placed in the f ir st location of the next to last set of KMAX 
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locations, and the coefficient associated with the +1 index 
should be placed in the first location of the last set of 
KMAX locations. 

To illustrate the effect of the selection of har- 
monics on the imperfections, equation (30) was solved using 
only the 13 selected harmonics. The resulting circumferen- 
tial imperfection rotation representation at station 10 
appears as curve B in Figure 18. To determine the effect of 
the inability to include all nine harmonics deemed important 
to the shell response with a fixed support, an alternate 
selection of harmonics was made for this case (see Table 
VII) . A complete SATANS-II analysis was performed with 
this alternate selection and the results were found to 
differ from those using the primary harmonic selection by 
approximately 2%. 

3 . Results 

Six analyses (boundary conditions A, B and C with 
clamped and simple supports) were made for each shell where 
the imperfections were reduced by a factor of 1000. These 
'minute' imperfections acted only as triggers for the modal 
deflections and the resulting solutions gave the shell's 
consistent critical bifurcation buckling loads. These 
results are presented in Table VIII under the category of 
minute imperfections. Selected cases were analyzed where 

the minute imperfection triggers were replaced by a 10 -7 psi 

« 

internal pressure trigger and the results were identical. 
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TABLE VIII 



BUCKLING LOADS, MINUTE VERSUS FULL IMPERFECTIONS 



SUPPORT 


DEGREE 

OF 

IMPERFEC- 

TION 


BOUNDARY 

CONDITION 


SHELL 


A- 14 


SHELL A-13 


BUCKLING 

LOAD 

(lb/in) 


REDUC- 

TION 

(%) 


BUCKLING 

LOAD 

(lb/in) 


REDUC- 

TION 

(%) 


C 




A 


47.20 




35.60 




L 


Minute 


B 


51.90 




50.60 




A 




C 


56.00 




55.10 




M 














P 














E 




A 


43.00 


8.90 


* 


- 


D 


Full 


B 


46.60 


10. 20 


40.70 


19.55 






C 


49.00 


12.50 


40.15 


22.70 






A 


48.00 




46.60 




S 


Minute 


B 


51.10 




50.60 




I 




C 


54.10 




55.30 




M 














P 














L 




A 


42.80 


10.8 


34.30 


26.30 


E 


Full 


B 


46.10 


9.8 


38.00 


24 .90 






C 


47.90 


11.5 


42.7 


22.80 



* insufficient data. 



The same six analyses for each shell were also made 
using the full imperfections. Results for the final load at 
which a converged solution could be obtained are also shown 
in Table VIII along with the percent reduction in the buckling 
load for each case. 

Since the full imperfection is introduced at the 
first load step as a pseudo load, the program requires con- 
siderably more iterations than normal to obtain the first 
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solution. The subsequent solution usually proceeds in a 
completely normal fashion. However, convergence difficul- 
ties may be encountered in cases where the imperfection 
level is high. This was the case with shell A-13, as some 
of its imperfection deviations exceeded 3 % shell wall thick- 
nesses, and fluctuated over 6 thicknesses in a 0 arc of only 
53°. For shell A-13, a 70% initial imperfection level was 
the maximum for which a converged solution could be obtained 
on the first load step. Consequently, buckling loads were 
computed with a maximum of 80%, 85%, 90% and 95% of the full 
imperfection data. For the 80% imperfection, 70% was intro- 
duced at load step one and 10% was added at load step two. 
For the three other solutions, 5% increments were added at 
various load steps. Since SATANS uses the total value of 
the variables and not the incremental values, each inter- 
mediate solution is correct for its given load and imper- 
fection level, but the load-deflection characteristics are 
invalid. An imperfection level of 95% was the absolute 
maximum for which a converged solution could be obtained for 

for shell A-13, and a short 5% extrapolation was made to ob- 

5 

tain the data presented in Table VIII for this shell. This 
procedure is illustrated for the shell with a simple support 
in Figure 19. A comparison of load steps and iterations for 
both shells under these conditions is presented in Table IX. 

^For shell A-13 with clamped boundary condition A, a con- 
verged solution was not obtained for an imperfection level 
above 85%. This condition is therefore not included in the 
results. A similar problem was encountered by Dantone /Ref. 
5/ with this boundary condition. 
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55 — 




30 — 



extrapolations to the ioo% 

IMPERFECTION LEVEL 



100 95 90 85 80 75 

IMPERFECTION LEVEL, % 

CRITICAL BUCKLING LOAD vs IMPERFECTION LEVEL 

FIGURE 19 
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TABLE IX 



LOAD STEPS VERSUS ITERATIONS 



SHELL A-14 


SHELL A-13 


LOAD 


INCREMENTAL 


ITERATIONS 


LOAD 


INCREMENTAL 


ITERATIONS 


STEP 


IMPERFECTION 


REQUIRED 


STEP 


IMPERFECTION 


REQUIRED 




ADDITION 






ADDITION 




1 


100% 


7 


1 


70% 


24 


2 




2 


2 


10% 


6 


3 




2 


3 


5% 


16 


4 




3 


4 




12 


5 




3 


5 




9 


6 




3 


6 




3 


7 




3 


7 




8 


8 




3 


8 




6 


9 




3 


9 


5% 


19 


10 




3 


10 




3 


11 




3 


11 


5% 


60 


12 




4 


12 




58 


13 




4 


13 




18 


14 




4 


14 




16 


15 




5 


15 




16 


16 




7 


16 




12 


17 




10 








18 




16 


1 st load step size reduction 


n 4* , 

1 load step size reduction 


17 




4 


19 




3 


2 nd ]_ oac j s tep size reduction 


20 




2 














18 




3 


2 n< ^ load step size reduction 


19 




2 


21 




3 


Trl 

3 load step size reduction 


r d 

3 load step size reduction 




termination 


22 




2 




(P cr =40.6 lb/in) 




termination 










(P cr =4 6 . 0 lb/in) 









Note: This data corresponds to boundary condition B, simple 

support . 
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D. DISCUSSION 



The SATANS-I I investigation of the buckling load of the 
perfect cylinder (discussed in section V.C.l) indicates that 
a numerical analysis using 35 stations results in a higher 
predicted buckling load than the correct analytical buckling 
load. Consequently, the SATANS-II results for the predic- 
tion of the critical buckling load of the imperfect cylinder 
are expected to be high. Table X compares the results of 
the SATANS-II analysis (column A) with the actual experi- 
mental results of Ref. 6 (column E) and shows this to be 
the case. The use of more stations to lower the predicted 
buckling load can be accomplished only by either decreasing 
the number of harmonics considered in the analysis or by 
increasing overall program size. Both of these choices are 
impractical at the present time since the number of harmon- 
ics already in the analysis appears to be low considering 
the number required for shell response and the imperfection 
data, and the present program core storage requirement is 
large. 

It is possible, however, to apply the SATANS-II results 
for the percent reduction in the buckling load due to the 
imperfections (see Table VIII and Column B, Table X) to a 
more accurate buckling load of the perfect shell from an 
analysis like that of Almroth to get a more accurate estimate 
of the buckling load of the imperfect shell. When the per- 
cent reduction from the consistent minute imperfection 
buckling load to the consistent full imperfection buckling 
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load, determined by SATANS-II, is applied to the accurate 
consistent bifurcation buckling load of Almroth (column C, 
Table X) , the result is a prediction of the actual buckling 
load of the imperfect shell. These imperfect shell predict- 
ed buckling loads are shown in column D of Table X and are 
seen to be very close to the actual buckling loads of the 
imperfect shells, given in Column E of Table X. 

The percent reductions in buckling loads determined 
by SATANS-II are considered to be applicable to the accurate 
consistent bifurcation buckling loads of Almroth since the 
imperfection characteristics are seen to be concentrated in 
the lower ten harmonics (see Table VI) , where shell response 
does not appear to be largely affected by the number of 
meridional stations in the analysis (see Figures 14 and 15) . 

E. PROPOSED ANALYSIS PROCEDURE 

When the SATANS-II predicted buckling load of the perfect 
shell is higher than the analytical prediction, as is the 
case with the extremely thin-walled cylinder, the following 
procedure is proposed. 

1. Determine the consistent bifurcation buckling load 
as close as possible. 

2. Determine the buckling load with SATANS-II using 
minute imperfections. 

3. Determine the buckling load with SATANS-II using 
full imperfections and compute the percent reduction from 
step 2. 

4. Apply this reduction to the buckling load of step 1. 
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This procedure is expected to work only when the result 



of step 2 is reasonably close to that of step 1. 
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VI. CONCLUSIONS 



A. SATANS-I 

The computer program SATANS-I has been constructed from 
SATANS. It is a more user-oriented version of the program 
by Ball, whose primary modification is the introduction of 
an equipment-independent graphical output ability. 

B. SATANS -I I 

SATANS -I I is a further modification of SATANS where the 
restriction that all applied loading and initial conditions 
be symmetrical about some meridional plane has been removed. 
In addition, SATANS-II possesses the capability of including 
experimentally-determined naturally-occurring imperfections 
or malformations induced during fabrication, mathematically 
described imperfections of any type, frequency or location, 
or graphically described imperfections due to damage suffer- 
ed after fabrication. The imperfections are introduced into 
SATANS-II by using the imperfection deviations along discrete 
circumferences to determine the associated meridional and 
circumferential imperfection rotations at the mesh points, 
performing a Fourier expansion on these rotations and then 
use of the Fourier coefficients. In cases where the impor- 
tant shell response harmonics and the harmonics required for 
imperfection representation exceed the limits of SATANS-II, 
a selection must be made. 
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Instructions for the use of both SATANS-I and SATANS-II 



are contained in Appendix C and complete listings of both 
programs are given in other appendices. SATANS-I requires 
210 k-bytes (53k words) of core on an IBM-360/67 computer 
using 'single precision,' while SATANS-II requires 426 
k-bytes (107k words) of core. Imperfect shell analyses 
performed by SATANS-II in this report used a 35 station 
meridian and 25 circumferential modes. Execution times 
were an average of 60 minutes, although execution times 
of between 15 and 30 minutes can be expected for analyses 
using 15 or less modes. 

C. IMPERFECT SHELLS 

SATANS-II has demonstrated an ability to perform a 
numerical structural analysis on the arbitrarily imperfect 
shell of revolution. Imperfections of over three shell wall 
thicknesses in depth have been handled. Increasing solution 
difficulty relates directly to increasing imperfection 
complexity and increasing shell instability. In those cases 
where extremely thin-walled shells result in perfect (imper- 
fection-free) buckling loads from SATANS that are higher 
than analytically determined buckling loads, a procedure has 
been devised where the percent reduction due to the intro- 
duction of full imperfections can be applied to the analyti- 
cally determined buckling load, and an imperfect shell 
buckling load predicted. These predictions have been shown 
to be very close to the actual imperfect shell buckling 
loads for the two specific axially loaded cylinders considered. 
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ADDENDUM 



A DIGITAL COMPUTER STUDY OF THE EFFECTS OF VARIOUS ARTIFICIAL 
BOUNDARY CONDITIONS AT THE POLE OF A SHALLOW SPHERICAL SHELL 

I. INTRODUCTION 

The design of many shell of revolution structures re- 
quires continuity of the structure at one or both poles. 

At such a pole, the radius becomes zero and special 
conditions must be imposed on the deflections and the stress 
resultants to assure finite strains and equilibrium. As a 
consequence of the difficulty of incorporating these 
conditions into a structural analysis computer program, 
many researchers model such a structure as having a small 
hole or a small rigid plug at the pole. This addendum 
investigates the applicability of the various artificial 
pole conditions by examining their effect on the stability 
of a clamped spherical cap subjected to axisymmetric and 
asymmetric pressure loads. 

Various polar artificialities are discussed and the 
boundary equations associated with their use are given, as 
are those for the continuous pole. A determination of the 
applicability of the various poles is made by comparisons 
of their effect on the shell buckling loads and by comparing 
the errors in the normal displacement and meridional mem- 
brane force distributions incurred by their use. 
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The computer program SATANS , developed by Ball /Ref. 37/ 
was used to conduct this study and selected results of the 
SATANS analyses are checked for validity with the computer 
programs BOSOR /Ref. 15/ and STAGS /Ref. 15/. SATANS has a 
special routine that incorporates the proper conditions to 
be enforced at a pole. BOSOR also has this special ability 

g 

but STAGS is not yet so equipped. 

SATANS, a one-dimensional finite difference equilibrium 
formulation based on Sanders' equations, was run on an 
IBM-360/67 computer under OS/MVT. BOSOR, a one-dimensional 
finite difference energy formulation based on the Donnell 
equations, was run on a UNIVAC-1108 computer under EXEC-VIII, 
while STAGS, a two-dimensional finite difference energy 
formulation was run on a CDC-6600 computer under COMPASS-12. 



6 



As of February, 



1972. 
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II. STATEMENT OF THE PROBLEM 



A. GEOMETRY AND LOADING 

The structure investigated is the clamped shallow 
spherical cap whose geometry can be specified by the single 
nondimensional parameter X, where 

X = 2(3(l-v 2 ))^(H/h) 33 (AD-1) 

and where H is the rise of the shell and h is its thickness. 
For this study, Poisson's ratio was taken as 0.3, Young's 
modulus as 30 x 10 6 psi and h as 1.0 inches. Shells of 
X = 6 and X = 18 were investigated however only the results 
of the X = 6 shell are given in detail. 

The classical buckling pressure of a complete sphere is 
given as q Q , where 

q Q = 2Eh 2 /R s /(3(l-v 2 )} J5 (AD-2) 

and where R g is the radius of curvature. For this study, 
two loading conditions were used. The cap was subjected to 
a uniform normal pressure (symmetrical load) and a uniform 
normal pressure over half the cap (asymmetrical load) . 

These loadings are shown in Figure 20. In each case, the 
magnitude of the applied pressure is q and the buckling 
results are presented as the ratio q/q Q . 
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X=6 SHALLOW SPHERICAL CAP 
GEOMETRY AND LOADING 

FIGURE 20 
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B. POLE CONDITIONS 



A shell is considered to have a pole if the structure is 
continuous through its axis of revolution. At a pole, the 
radius is zero and special conditions must be imposed on 
deflections u, v and w, and on the meridional moment m g to 
assure finite strains and equilibrium. The computer program 
SATANS is equipped with the capability to impose these 
special boundary conditions. A full derivation of the 
required conditions is given in Appendix D of Ref. 3. The 
results of that derivation are presented below. 

For the zeroeth Fourier harmonic (n=0) : 



W " = m s - = o 

rn a = m 
so s 

t Q = t s = b(i-v) (V+w s w+^e s )-t T 



(AD-3) 



For the first Fourier harmonic (n=l) 



( u ± v ) = W = u" = 0 



m s m 0 ra S 0 ~ 0 

S = S - So = 0 



(AD-4 ) 



When a double sign (±) appears, the upper sign applies for 
an initial pole and the lower sign applies for a final pole. 
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For the second harmonic (n=2) : 



u = v = w = m^ = 0 
s 

fc 0 = _t s = ±t s0 
m 0 = ~ m s = ±m s 0 
t s = b(l-v) (u'+izBg) 



(AD-5) 



For all higher harmonics (n>2) : 
u = v = w = 0 

m s = m 0 = ra s 0 = 0 (AD “ 6) 




For rotations at the pole, for the first harmonic (n=l) 



<P S = -w'+w s u 

<P Q = :*4> s 

<P = %[2v'±u') 

and for all other harmonics, 

<{>3 = 4> 0 = 0 

<f> =h(2v'±nu') 



(AD-7) 



(AD-8) 



Numerous artificial boundary conditions can be substi- 
tuted for the correct conditions at the pole. For example, 
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a small hole can be placed at the pole and can exist with 
the boundary completely free, or have conditions imposed 
such that the meridional rotation <f s and deflection u are 
restrained. A meridional line load N g of magnitude suffi- 
cient to replace that 'lost' due to the hole can be applied 
Another artificial pole condition is that best described as 
a small rigid plug. Here, both the meridional as well as 
the circumferential deflection is restrained along with the 
meridional rotation. A normal line load -P whose total 
equals that loading 'lost' due to the hole or plug can be 
applied at the boundary of the hole or plug. 

Preliminary numerical studies showed that the line load 
ings effected a noticeable change in the solution only in 
the immediate area of the boundary, hence, the line loads 
were not considered further. Preliminary studies also 
showed that the cap behavior exhibited only minor changes 
when the small hole with a partially restrained edge was 
substituted for that with a totally free edge. This 
investigation, therefore, concentrated on two artificial 
boundary conditions at a pole and on the continuous pole. 
The artificial pole conditions used and their associated 
mathematical expressions are: 

1. The totally free hole: 

N s. = 0 ; N s0 = 0 ; Q s = 0 ' M s = 0 (AD-9) 

2. The rigid plug: 

U=0;V=0;Q S = 0;4> S =0 (AD-10) 
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c. TYPES OF ANALYSIS 

Three types of buckling analyses were conducted on the 
cap under each of the three pole conditions and these were: 

1. The nonlinear axisymmetr ic analysis of the symmet- 
rically loaded clamped spherical cap. 

2. The bifurcation buckling analysis of the symmetri- 
cally loaded clamped spherical cap with nonlinear prebuckling 
axisymmetric displacements. 

3. The nonlinear collapse analysis of the asymmetrically 
loaded clamped spherical cap. 

In the artificial pole cases, hole and plug radii were 
taken as 0.17%, 1.70% and 3.30% of the shell meridian. The 
SATANS analyses used a 100 station meridian for the axisym- 
metric cases, a 66 station meridian with three harmonics for 
the bifurcation buckling cases and a 40 station meridian 
with five harmonics for the nonsymmetric collapse cases. 
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III. DISCUSSION 



The results of the SATANS analyses are presented in 
Table XI for all cases. Also presented are the normal 
deflections of a point located at 47.5% of the meridian 
which is near the maximum displacement for each case. 

A. THE X = 6 SHALLOW CAP 

1 . The Ax i symmetric Analysis 

In Table XI, the buckling load is seen to be basi- 
cally unaffected by the altered pole conditions. The normal 
displacement, however, is seen to show a marked sensitivity 
to the pole conditions. Figure 21 shows this sensitivity. 

In the figure, the upper curve is the deflection shape of 
the cap at incipient buckling with a continuous pole, and 
the lower curves represent the percent error induced in 
the normal displacement at incipient buckling by the intro- 
duction of artificial poles. The 0.17% plug and the 0.17% 
hole behaved identically. Holes, in general, produced an 
average 20% error in the normal displacement and this error 
existed over the entire shell . The plugs, preserving the 
continuity of the structure, produced an average 10% error, 
but again, this error extended over the entire shell. 

Figure 22 shows the behavior of the meridional force N s for 
the continuous pole case, and the N s error due to the 
artificial pole conditions, also at incipient buckling. 

Large errors existed in all cases over the first third of 
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TABLE XI 



SHALLOW CAP (A=6) BUCKLING LOADS AND DISPLACEMENTS AT A POINT 



ANALYSIS 

TYPE 


POLE 

CONDITION 


MERIDIONAL 

STATIONS/ 

HARMONICS 


q/q Q 


NORMAL 
DEFLECTION 
AT 47.5% 
MERIDIAN 
( inches) 


Ax i symmetric 
Buckling 

Symmetric 

Load 


Pole 


3 0/1 
100/1 


0.9950* 

1.0610 

1.0200 


1.3720 

1.3200 


-0.17% 
Hole-1. 70% 
-3.30% 


100/1 

100/1 

100/1 


0.9720 

0.9690 

0.9620 


1.1600 

1.0560 

1.0290 


-0.17% 
Plug-1. 70% 
-3.30% 


100/1 

100/1 

100/1 


0.9720 

0.9770 

0.9900 


1.2570 

1.1630 

1.2290 


Bifurcation 

Buckling 

Symmetric 

Load 


Pole 


3 0/3 
66/3 


0.7750* 

0.7851 

0.7752 


0.6386 

0.6188 


-0.17% 
Hole-1. 70% 
-3.30% 


66/3 

66/3 

66/3 


0.7719 

0.7835 

0.7392 


0.6356 

0.6562 

0.6264 


-0.17% 
Plug-1. 70% 
-3.30% 


66/3 

66/3 

66/3 


0.7719 

0.7726 

0.7713 


0.6356 

0.6438 

0.6513 


Collapse 

Analysis 

Asymmetric 

Load 


Pole 


4 0/5 


0. 7100+ 
0.5391 


0.9830 


-0.17% 
Hole-1. 70% 
-3 .30% 


4 0/5 
4 0/5 
40/5 


0.5363 

0.5389 

0.5962 


1.0780 

1.0060 

1.4110 


-0.17% 
Plug-1. 70% 
-3.30% 


40/5 
4 0/5 
40/5 


0.5978 

0.6290 

0.7576 


0.9478 

0.8100 

0.7248 



* Huang /Ref. 16/ 
t Famili and Archer /Ref. 17/ 



99 




700 
saq'oui 1 



lN3tM33VldGia IVVMdON 




0 J 3<j 



!N3lM3DV~ldSia IVWdON 



CL 

< 

o 



CD 

it 

CD 

CD 

>- 

-J 

< 

2 

< 

a 

cl 

\r~ 

UJ 

2 

>- 

CD 

X 

< 

cl 

o 

cl 

CL 

UJ 

H 



cvj 

UJ 

CL 

ZD 

CD 

U_ 



UJ 

2 

UJ 

o 

< 

-J 

CL 

CD 

Q 

-J 

< 

2 

CL 

O 



100 




0L 

< 

O 

CD 

ii 

CO 

CO 

>- 

-J 

< 

z 

< 

o 

a: 

h* 

LlJ 

CO 

X 

< 

of 

o 

rr 

o: 

uj 

LU 

o 

cc 

o 

n. 

_i 

< 

z 

o 

Q 

OC 

UJ 

2 



CM 

04 

LU 

CL 

3 

O 

u_ 



101 



the meridian and were an average of 10% or less for the 
remainder of the meridian. 

2. Bifurcation Buckling Analysis 

In this analysis, the shell buckling load appears to 
be insensitive to artificial pole conditions for the case of 
the plugs, and shows some sensitivity in the case of the 
holes. This is surprising since bifurcation buckling of the 
X = 6 cap occurs in the Cos (20) mode. While equations (AD-8) 
show that no rotations exist at the pole in this mode, rota- 
tions certainly do exist at 1.7% and 3.3% meridian due to 
the buckling mode. The equations for the rigid plug's 
boundary conditions (equations (AD-10) ) specifically pre- 

vent this rotation, while those for the hole (equations 
(AD-9) ) do not. 

Examination of Figures 23 and 24 show that the dis- 
placement and force distributions of the plugs at incipient 
buckling to be less erroneous, again, than the poles. In 
Figure 23, an average 10% displacement error occurs over the 
entire shell for the plugs, while holes of 1.7% and 3.3% 
meridian cause large errors over much of the shell. This 
trend is repeated in Figure 24, although the errors are not 
as large. From this analysis, the indications are that the 
1.7% hole is a poor artificial pole and a 3.3% hole is 
unacceptable as it causes a 4.6% error in buckling load and 
over 40% error in some regions in both normal displacement 
and meridional force distributions. Larger hole sizes would 
apparently increase these errors. 
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NORMAL DISPLACEMENT ERROR, BIFURCATION BUCKLING ANALYSIS, A= 6 CAP 

FIGURE 23 
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MERIDIONAL FORCE ERROR, BIFURCATION BUCKLING ANALYSIS, X = 6 CAP 

FIGURE 24 . 



3 . Non symmetric Collapse Analysis 

The results for the collapse analysis are also pre- 
sented in Table XI, and in Figures 25 and 26 for the 0=0° 
meridian, which is symmetrically located under the loaded 
region. In Table XI, it is shown that the buckling loads of 
the cap with rigid plugs are quite erroneous (10.8%, 16.6% 
and 41.4% for the 0.17%, 1.7% and 3.3% plugs, respectively) 
while those for the small holes are quite reasonable. In 
Figures 25 and 26, it is seen that the small holes cause an 
average error of 10% in both the normal displacement and 
meridional force distributions at incipient buckling over 
the entire shell, while the errors caused by the rigid 
plugs is large. 

In an asymmetrically loaded shell, large rotations 
will exist at the pole and the restrictions imposed by the 
rigid plug appear to make it an unacceptable artificial pole 
for use in this case. 

Figures 25 and 26 also indicate that a 1.7% hole is 
the best artificial pole for this case. Assuming such a 
'best' hole size exists for another given shell or loading 
condition, an investigation similar to that presented here 
would have to be performed to find it. This is an inherent- 
ly expensive and time consuming process and there is no 
guarantee that there is a 'best' hole size without its 
radius becoming so small that mathematical difficulties 
arise during the solution process. Consequently, the use' of 
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NORMAL DISPLACEMENT ERROR, NONSYMMETRIC COLLAPSE ANALYSIS, A= 6 CAP 

FIGURE 25 
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artificial pole conditions in nonsynunetr ic problems must be 
approached with caution. 

B. THE X = 18 SHALLOW CAP 

A deeper shell with X = 18 was investigated in a similar 
fashion and produced similar results. The shell buckling 
load was negligibly affected (except in those cases already 
mentioned for the X = 6 shell) but displacement errors of up 
to 30% and meridional membrane force errors of an average 
10-15% were noted. In general, the X = 18 shell appeared 
less sensitive to the type of artificial pole (except in the 
asymmetric case) , but as sensitive to any artificial pole as 
the shallower X = 6 shell. 

C. COMPARISONS WITH OTHER RESEARCHERS 

As indicated in the introduction, various SATANS results 
were checked for validity with BOSOR and STAGS. These com- 
parisons are presented with others in Table XII and very 
favorable agreement is seen to exist. 
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SHALLOW CAP (X=6) BUCKLING LOADS COMPARISON 




* 
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This result and the following three were taken from a discussion of the effects of the 
number of stations and harmonics on the buckling loads of the clamped spherical cap, 
contained in Ref. 18. 



IV. CONCLUSIONS 



The stability of the sy^etrically loaded clamped shal 

low spherical shell appears to be fa,'rlv • 
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5. In general, the greater the displacement, or the 
higher the membrane force that would exist at a pole, the 
greater the error introduced by the use of an artificial 
pole . 
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APPENDIX A 



NONLINEAR TERMS 

The nonlinear terms of equations (15) are given below. 
Terms involving the initial imperfection rotation terms 
along have been omitted as they cause no straining. 
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APPENDIX B 



EIGENVALUE PROBLEM MATRIX ELEMENT DEFINITIONS 

The elements of the /A_/ and /B_/ matrices used in the 
eigenvalue solution are defined below. 



i i 1 



= G 



1 1 



- e 11 q : 



^12 Ej 2 Q ' 



‘1 3 



= F 13 Q 



A 2 x 
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A 23 = 
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E 23 Q 2 ~ 
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and where all elements of the 4x4 matrices / E_7, / F / 
and / — G 7 are zero except for the following, defined as 



En = b 

E 22 = 3s(l-v)b + 9 ( 1-v) d/ ( 8 Rq ) 

E 23 = 3 (1-v) dn/ (2R 2 ) 

E 32 = E 23 (B-3) 

E 3 3 = 2 (1-v) dn 2 /R 2 
E 34 = 1 

E 4 3 = ~d 

F 1 2 = h (1-v) bn/R 0 - 3 (l-v)dn/( 8 R 3 ) 

F 1 3 = bv/R A - Ml-v)n 2 /Rn 

6 6 (B-4) 

F 2 1 = -F 1 2 
F 3 1 = “F 1 3 

G 11 =~h (1-v) bn 2 /R 2 - (1-v) dn 2 / ( 8 Rg) 

G22 =-bn 2 /R 2 - (1-v) (1+v) dn 2 /Rg 
G 23 = -bn/Rg - (1-v) (l+v)dn 3 /RQ 
G 2 4 = -vn/R 2 
G 3 2 = G 2 3 
G33 = -vn 2 /R 2 
G 42 = vdn/R 2 
G 4 3 = vdn 2 /R 2 
G 4 4 = —1 
and where 

b = { ( 1 -v) ( 1 +v ) } _1 
d = (b/ 12 ) (h/r ) 2 



(B-5) 



(B-6) 
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APPENDIX C 



INPUT DATA GUIDE FOR SATANS-I AND SATANS-II 



Card 


Columns 


Format 


Item 


Example 


Meaning 


1 


1-72 


18A4 


TITLE 


- 


Enter any 72 characters 


2 


1- 5 


15 


NO 


1 


The problem number, 
0<NO<10000. 


2 


6-10 


L5 


$DYNAMC 


F 


For a static analysis, 
set $DYNAMC = F. 










T 


For a dynamic analysis, 
set $DYNAMC = T. 


2 


11-15 


15 


IMODE 


0 


For no modal output data. 










1 


For modal output data. 


2 


16-20 


15 


NDIMEN 


0 


Dimensional output data. 










1 


Nondimensional output. 


2 


21-25 


15 


NTHMAX 


8 


Summed solution will be 
printed at NTHMAX merid- 
ians, 0<=NTHMAX<=36 . 


2 


26-30 


15 


IFREQ 


3 


Solution will be printed 
at the first station, 
every subsequent IFREQ 
station and the last sta- 
tion, 0<IFREQ<=KMAX . 


2 


31-35 


15 


IPRINT 


2 


Every IPRINT converged 
solution will be print- 
ed . 


2 


36-40 


15 


IBCINL 


-1 


If the shell has a pole 
at the first station. 










0 


If the shell has no pole 
at the first station. 


2 


41-45 


15 


IBCFNL 


-1 


If the shell has a pole 



at the last station. 

If the shell has no pole 
at the last station. 
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Card Columns Format Item 



Example Meaning 



2 46-50 15 KMAX 35 

2 51-55 15 MNMAX 7 



2 61-65 15 LSMAX 1 

99 

750 

2 66-70 15 LCHMAX 3 

0 

2 71-75 15 ITRMAX 1 

30 

2 76-80 15 IC 0 

1 



Number of meridional sta- 
tions . 

Note: KMAX<201 for SATANS 
-I without plots and 
KMAX<101 for SATANS-I 
with plots or for SATANS 

-II. 

Number of series coeffi- 
cients used to describe 
the initial conditions, 
pressure and thermal 
loads (and initial imper- 
fections if using SATANS 
-II) . 

Note: For SATANS-I: 

0<MAXM<11 and 
KMAX*MAXM< 2 01 . 
For SATANS-I I: 
0<MAXM< 2 6 and 
KMAX*MAXM< 951. 

For a linear analysis. 

Use many load steps for 
a nonlinear static analy- 
sis . 

For a dynamic analysis, 
LSMAX is the number of 
time increments, where 
LSMAX=T max /AT. 

The number of load step 
size reductions in a 
static analysis, recom- 
mended range = 2-4. 

For a dymanic analysis. 

For a linear analysis. 

The number of iterations 
at a load or time step. 
For a nonlinear analysis, 
suggested range = 10-30, 
up to 50 for special 
cases . 

Initial conditions. Set 
to 0 for a static analy- 
sis, or for a dynamic 
analysis where the shell 
is at rest at T=0. 

For a dynamic analysis 
with initial conditions. 
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Card 


Columns 


Format 


Item 


Example 


Meaning 


3 


1-12 


E12.3 


NU 


0.3 


Poisson's ratio, v. 


3 


12-24 


El 2 . 3 


SIGO 


1000.0 


Reference stress level, 
Oq • 










1.0 


If the input data is 
dimensional . 


3 


24-36 


E12.3 


ELAST 


. 3E8 


Reference modulus of 
elasticity, E 0 . 










1.0 


If the input data is 
dimensional . 


3 


u> 

I 

00 


E12.3 


TKN 


. 4E-2 


Reference thickness, h Q 










1.0 


If the input data is 












dimensional . 


3 


49-60 


E12.3 


CHAR 


8.16 


Characteristic shell 
dimension, a. 










1.0 


If the input data is 
dimensional . 


3 


61-72 


E12.3 


TEEO 


0.0 


If a static analysis. 










. 996E-5 


Reference time, T 0 . 


4 


1-12 


E12.3 


DELOAD 


0.2 


For a static analysis. 



DELOAD is the load incre- 
ment. It remains un- 
changed until the solu- 
tion fails to converge 
in ITRMAX iterations, 
when it is reduced by a 
factor of five. A maxi- 
mum of LC1IMAX such re- 
ductions will occur. 

. 18 23E-6 For a dynamic analysis, 
DELOAD is the nondimen- 
sional time increment. 

4 13-24 E12.3 EPS 0.01 The convergence criterion 

recommended range of 
0. 0KEPSC0.001. 



Card 4A is only included for a SATANS-II analysis. 

4a 1- 5 15 JUMP I For an analysis using 

single series expansions. 
For an analysis using 
double series expansions. 



2 



Card Column Format Item 



Example Meaning 



4A 5-10 15 MPERFS 0 An analysis without im- 

perfections . 

1 An analysis with imperfec- 
tions. 

Note: if JUMP=2 , MPERFS 
may be 0 or 1 . 

if JUMP=1, MPERFS 
must be 0 

if MPERFS=1, JUMP 
must be 2. 



Include as many cards 5 as necessary to specify NTHMAX 
meridians. If NTHMAX equals 0, omit card 5. 

5 1-72 6E12.3 10.0 A list of circumferential 

coordinates 6, in degrees 
and tenths, where the 
solution printout is de- 
sired. The list must have 
NTHMAX entries. 



If IBCINL=-1, omit cards 6 through 14. If IBCFNL=-1, omit 
cards 15 through 23. Cards 6 through 23 describe the 
boundary conditions at the first, and then at the last sta- 
tion. The boundary conditions exist on the total variables, 
not on the individual harmonics. Loadings applied through 
specification of boundary conditions are taken in the zeroeth 
harmonic (n=0) only, as the column matrix {£} is set to zero 
for harmonics greater than zero. The boundary conditions are 
dimensional. The format of cards 6 through 23 is 4E16.8. 



Card 6,15 Card 



0(1,1) 


0(1 


0(2,1) 


0(2 


0(3,1) 


0(3 


0(4,1) 


0(4 



7,16 Card 8,17 
2) o(i,3) 

2) 0(2,3) 

2) 0(3,3) 

2) 0(4,3) 



Card 9, 18 



0(1,4) 




f \ 

N s 


0(2,4) 




N s e 


0(3,4) 




Q s 


0(4,4) 




*s 


1 




L 



+ 
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Card 10/19 Card 11,20 Card 12,21 Card 13,22 



Card 14,23 



r 

A (1 , 1) 


A (1 , 2) 


A Cl, 3) 


A(l, 4) 




u 




' 

A(l) 


A (2,1) 


A (2, 2) 


A (2 , 3 ) 


A (2,4) 




V 




M2) 


A (3,1) 


A (3 , 2) 


A (3 , 3 ) 


A (3 , 4 ) 




w 




M3) 


A (4 , 1) 


A (4 , 2) 


A (4 , 3 ) 


A (4 ,4) 




M 

s 




M4) 


k 






4 




k J 




V J 



Card 24 is: 

1. Included for a SATANS-I static analysis. 

2. Included but blank for a SATANS-I dynamic analysis. 

3. Omitted for a SATANS-II analysis. 



Card Column Format Item Example Meaning 



24 1- 2 L2 $PLOTS F Indicates plots are NOT 

desired . 

T Indicates plots are 

desired . 



24 3- 4 L2 $MODAL F Indicates plots are for 

summed solutions only 
T Indicates plots are for 

modal solutions only. 



For the remainder of card 24 entries, 0 indicates that no 
plots are desired for the particular item, and 1 indicates 



that they are desired. All graphs are plotted as the in- 
dicated item versus the station number. If a complete plot 
is desired, insure IFREQ=1. 


24 


5- 6 


12 


I RADI I 


1 


Plot the radii as computed 
by subroutine GEOM. 


24 


7- 8 


12 


I GAMMA 


1 


Plot p'/p as computed by 
subroutine GEOM. 


24 


9-10 


12 


IOMEGS 


1 


Plot a) as computed by 
subroutine GEOM. 


24 


11-12 


12 


IOMEGT 


1 


Plot u) fl as computed by 
subroutine GEOM. 


24 


13-14 


12 


IDEOMS 


1 


Plot ojg as computed by 
subroutine GEOM. 


24 


15-16 


12 


IBSTIF 


1 


Plot the stiffness d as 



computed by subroutine 
BDB . 
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Card 


Column 


Format 


Item 


Example 


Meaning 


24 


17-18 


12 


IDSTIF 


1 


Plot the stiffness d as 
computed by the subroutine 
BDB . 


24 


19-20 


12 


IBBSTF 


1 


Plot the stiffness db/d? 
as computed by subroutine 
BDB. 


24 


21-22 


12 


IDDSTF 


1 


Plot the stiffness dd/d? 
as computed by subroutine 
BDB. 


24 


23-24 


12 


IPR 


1 


Plot the normal component 
of the pressure load. 


24 


25-26 


12 


IPS 


1 


Plot the meridional com- 
ponent of the pressure 
load . 


24 


27-28 


12 


IPT 


1 


Plot the circumferential 
component of the pressure 
load. 


24 


29-30 


12 


ITT 


1 


Plot the thermal load. 


24 


31-32 


12 


IMT 


1 


Plot the thermal moment. 


24 


33-34 


12 


IDTT 


1 


Plot d/d? of the thermal 
load . 


24 


35-36 


12 


I DMT 


1 


Plot d/d? of the thermal 
moment . 


24 


37-38 


12 


INS 


1 


Plot the meridional mem- 
brane force distribution. 


24 


39-40 


12 


INTH 


1 


Plot the circumferential 
membrane force distribu- 
tion 


24 


41-42 


12 


INSTH 


1 


Plot the meridio- 
circumf erential membrane 
force distribution. 


24 


43-44 


12 


IQS 


1 


Plot the transverse force 
distribution . 


24 


45-46 


12 


IMS 


1 


Plot the meridional moment 
distribution . 


24 


47-48 


12 


IMTH 


1 


Plot the circumferential 



moment distribution. 
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Card 


Column 


Format 


Item 


Example 


Meaning 


24 


49-50 


12 


IMSTH 


1 


Plot the meridio- 
circumf erential moment 
distribution . 


24 


51-52 


12 


IU 


1 


Plot the meridional dis- 
placement distribution. 


24 


53-54 


12 


IV 


1 


Plot the circumferential 
displacement distribution 


24 


55-56 


12 


IW 


1 


Plot the normal displace- 
ment distribution. 


24 


57-58 


12 


IPHIS 


1 


Plot the meridional rota- 
tion distribution. 


24 


59-60 


12 


IPHIT 


1 


Plot the circumferential 
rotation distribution. 


24 


61-62 


12 


IPHI 


1 


Plot the meridio- 
circumf erential rotation 
distribution. 



Insert imperfection data here for a SATANS-II analysis with 
imperfections. Insure format of the imperfection data is 
compatible with that specified in the user-written sub- 
routine IMPERF. 



25 1- 2 12 IRNAGN 0 Indicates this is the 

only run. 

1 Indicates another run is 

to be made. Add another 
complete set of data 
cards after this card is 
IRNAGN=1 . 
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THIS LISTING COMPRISES THE STORABLE PACKAGE OF SATANS-I AND CONTAINS THE 
FOLLOWING SUBROUTINES: 
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